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Introduction 


In  this  unit  some  models  are  discussed  that  attempt  to  describe  the  way  in  which 
populations  grow.  Such  models  are  used  in  demography — the  study  of  human 
populations— and  in  population  ecology— the  study  of  how  more  general 
biological  populations  interact  with  their  environment  and  with  other  populations. 

A  number  of  stochastic  models  for  the  growth  of  a  population  have  been  discussed 
in  earlier  units;  for  example,  branching  processes  were  studied  in  Unit  4  and  birth 
and  death  processes  in  Units  7  and  8.  Corresponding  deterministic  models  have 
also  been  discussed  briefly.  Most  of  the  models  in  this  unit  are  deterministic  rather 
than  stochastic.  The  first  model  discussed  in  Section  1  is  the  deterministic  analogue 
of  the  simple  birth-death  process,  and  is  called  the  exponential  model.  It  provides 
a  mathematical  expression  of  Malthus’  assertion  that:  ‘Population,  when 
unchecked,  increases  in  a  geometrical  ratio.’  The  second  model  incorporates  the 
idea  that  in  reality  there  is  a  limit  to  how  large  a  population  can  grow.  This 
model— the  logistic  model— may  be  regarded  as  a  modification  of  the  exponential 
model. 

In  both  these  models  the  rate  of  growth  of  the  population  is  assumed  to  depend 
only  on  the  current  size  of  the  population,  and  the  more  detailed  structure  of  the 
population  is  not  considered.  However,  more  accurate  descriptions  and  predictions 
of  population  growth  may  usually  be  obtained  if  factors  other  than  just  the 
population  size  are  taken  into  account.  Of  particular  importance  is  the  age- 
structure  of  the  population  that  is,  the  distribution  of  the  ages  of  the  members  of 
the  population— and  the  proportions  of  males  and  females  in  the  population:  two 
populations  of  the  same  size,  even  if  identical  in  other  respects,  may  have  very 
different  growth  rates  if  their  age-structures  and  sex-ratios  are  different.  To  take  an 
extreme  case,  if  a  population  consisted  predominantly  of  the  elderly,  then  the  birth 
rate  would  be  small  and  the  death  rate  large.  On  the  other  hand,  if  a  large 
proportion  of  the  females  in  a  population  were  in  the  fertile  age  groups,  then  we 
should  expect  the  birth  rate  to  be  relatively  large. 

Age-structure  may  be  taken  into  account  in  models  for  various  biological 
populations,  but  the  models  in  Sections  2  to  4  are  developed  with  human 
populations  in  mind.  There  is  an  abundance  of  data  concerning  the  age-structure 
and  birth  and  death  rates  of  human  populations:  records  of  births  and  deaths  have 
been  kept  in  many  countries  for  many  years,  and  information  on  the  size  and  age- 
structure  of  human  populations  is  readily  available.  Section  2  is  devoted  to  the 
study  of  some  of  the  fundamental  concepts  and  results  of  mathematical 
demography.  Data  on  the  age-structure  of  a  population  is  presented  in  the  form  of 
a  life  table;  this  gives,  for  various  ages  x,  the  proportion  of  individuals  in  a 
population  who  have  survived  to  age  x.  Life  tables  are  discussed  and  a  population 
model  based  on  the  life  table  is  developed;  this  is  used  to  investigate  the  way  in 
which  death  rates  depend  on  age. 

The  implications  of  this  model  for  the  growth  and  age-structure  of  a  population 
are  investigated  in  Sections  3  and  4.  In  Section  3,  a  population  whose  size  and  age- 
structure  do  not  change  with  time — a  stationary  population — is  studied;  and  in 
Section  4,  the  properties  of  a  population  whose  size  may  change  but  whose  age- 
structure  remains  constant — a  stable  population — are  discussed.  Finally,  in 
Section  5  the  age-structures  of  two  specific  human  populations  are  discussed 
qualitatively. 

Section  1  of  this  unit  draws  on  several  ideas  which  have  been  introduced  in  earlier 
units;  those  parts  of  Units  7  and  8  which  deal  with  obtaining  the  deterministic 
model  corresponding  to  a  birth  and  death  process  are  particularly  relevant.  The 
remaining  sections  do  not  depend  to  any  great  extent  on  earlier  parts  of  the 
course.  Section  2  contains  more  questions  than  do  the  other  sections;  you  will 
probably  find  that  it  will  require  most  study  time.  On  the  other  hand,  Section  5  is 
very  short  and  does  not  contain  any  questions,  so  you  may  wish  to  study  Sections 
4  and  5  in  one  session.  There  are  no  audio  or  video  components  associated  with 
this  unit. 


T.  R.  Malthus,  An  Essay  on  the 
Principle  of  Population  (Penguin, 
1970).  This  classic  work  on 
demography  was  first  published  in 
1798,  and  you  may  find  it  interesting 
reading. 
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1  Two  deterministic  models  for  the  growth  of  a 
population 


In  this  section  two  simple  models  for  the  growth  of  a  population  are  discussed  and 
their  properties  are  investigated.  The  first  of  these,  the  exponential  model,  you 
have  met  previously  as  the  deterministic  analogue  of  the  simple  birth-death 
process.  The  second,  the  logistic  model,  may  be  regarded  as  a  modification  of  the 
first;  in  this  model  the  growth  rate  at  any  time  depends  on  the  current  size  of  the 
population. 


The  exponential  model 

The  simple  birth-death  process  was  discussed  in  detail  in  Sections  3  and  4  of 
Unit  8.  For  this  stochastic  process,  the  overall  birth  rate  px,  and  the  overall  death 
rate,  vx,  are  proportional  to  the  size  of  the  population,  x;  that  is, 

Px  =  fix,  x  =  0,  1,  ...  , 

Vx  =  vx,  x=l,  2, ..., 
where  p  and  v  are  positive  constants. 


If  x  —  x(t)  is  the  size  of  the  population  at  time  t  in  the  corresponding  deterministic 
model,  then  the  differential  equation  for  this  model  is 
dx 

=  Px  —  vx,  t  >  0, 


and,  letting  0  =  p  —  v,  this  may  be  written  as 


(1.1) 


When  the  population  size  is  x,  the  rate  at  which  the  population  grows  is  Ox,  so  the 
1  dx 

per  capita  growth  rate,  -— ,  is  0.  The  parameter  0  is  called  the  intrinsic  per  capita 

growth  rate  of  the  population;  it  is  also  sometimes  referred  to  as  the  Malthusian 
parameter. 


If  the  population  size  is  x0  at  time  0,  then  Equation  (1.1)  has  the  solution 

x(f)  =  x0e°‘,  t>  0.  (1.2) 

If  0  >  0  (that  is,  if  j 6  >  v),  then  the  model  predicts  that  the  population  size  will 
grow  exponentially,  and  for  this  reason  it  is  often  referred  to  as  the  exponential 
model.  If  0  <  0,  then  the  model  predicts  that  the  population  size  will  decrease 
exponentially,  x(£)  -►  0  as  t  ->  oo,  so  the  population  decreases  in  size  but  never  dies 
out  completely.  And  if  0  =  0,  then  x(t)  =  x0  for  all  t  >  0. 

Writing  e9  =  </>  in  Expression  (1.2),  we  have 
x(t)  =  xo0',  t  >  0. 

So  after  each  unit  of  time,  the  size  of  the  population  is  multiplied  by  the  constant 
factor  (f>\  x(0)  =  x0,  x(l)  =  x0(J>,  x(2)  =  x0<f>2,  ...,  and  in  the  case  0  >  0,  that  is 
(f>  >  1,  this  shows  explicitly  the  equivalence  of  the  exponential  model  to  Malthus’ 
assertion  that  ‘population,  when  unchecked,  increases  in  a  geometrical  ratio’.  Note, 
however,  Malthus’  important  qualification  that  when  unchecked,  population 
increases  geometrically.  If  a  species  is  introduced  into  a  large  (but  finite)  habitat, 
then  initially  it  may  well  grow  exponentially,  but  eventually  the  limitations  in  food 
and  other  resources  imposed  by  the  habitat  will  cause  the  rate  of  growth  to 
decline.  To  put  it  another  way,  Expression  (1.2)  implies  that  if  0  >  0,  then  the 
population  grows  without  bound,  which  cannot  be  strictly  true,  since  no  finite 
habitat  can  support  an  unlimited  size  of  population.  We  may  suppose  that  any 
given  habitat  has  a  maximum  size  that  it  can  support;  this  is  called  the  carrying 
capacity  of  the  habitat  for  the  population. 


These  two  deterministic  models  are 
discussed  in  MST204  Unit  3. 


See  Unit  8,  Question  3.10. 


In  Question  3.10  of  Unit  8,  the 
behaviour  of  the  simple  birth-death 
process  was  compared  with  that  of 
the  analogous  deterministic  model. 
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How  can  the  exponential  model  be  modified  to  take  account  of  the  carrying 
capacity  of  a  habitat?  Do  any  of  the  models  that  have  been  discussed  in  earlier 
units  have  the  property  that  the  population  size  increases  with  time  but  can  never 
exceed  some  maximum  value?  In  fact,  one  of  the  models  that  you  met  in  Unit  10 
has  this  property;  this  is  the  deterministic  model  for  the  simple  epidemic  that  was 
discussed  in  Subsection  2.2  of  that  unit.  In  this  model,  y,  the  number  of  infected 
individuals  at  time  t  in  a  community  of  size  n  +  1,  is  the  solution  of  the  differential 
equation 

dy  _  (]y(n  +  1  -  y) 

dt  n 

The  number  of  infectives  increases  with  time,  becoming  closer  and  closer  to  n  +  1, 
although  the  whole  population  is  infected  only  after  an  infinite  time. 


The  logistic  model 


The  simplest  model  for  the  change  in  size  of  a  population  which  incorporates  the 
idea  of  a  carrying  capacity  is  equivalent  in  form  to  the  deterministic  model  for  the 
simple  epidemic.  It  is  expressed  in  the  differential  equation 


t  >0, 


(1.3) 


where  9,  the  intrinsic  per  capita  growth  rate,  and  the  parameter  k,  are  both 
positive;  this  model  is  called  the  logistic  model.  The  per  capita  rate  of  growth  when 
the  population  size  is  x  is  given  by 


a  decreasing  linear  function  of  x,  and  is  zero  for  a  population  of  size  k.  So  as  the 
population  changes  in  size,  the  rate  at  which  individuals  give  rise  to  new 
individuals  decreases,  approaching  zero  as  the  population  size  approaches  k.  If  a 
model  is  sought  for  a  population  which  is  growing,  it  will  also  be  assumed  that 

dx 

0  <  x  <  k,  so  that  the  rate  of  growth,  — ,  is  positive. 


Question  1.1 

(i)  For  0  <  x  <  k,  show  by  the  method  of  separation  of  variables  that  x,  the 
general  solution  of  the  differential  equation  (1.3),  satisfies 


log 


k-x 


=  —Ot  +  constant. 


GO 


If  the  population  size  is  x0  at  time  0,  where  0  <  x0  <  k,  show  that  the 
population  size  at  time  t  >  0  is  given  by 


x(f)  = 


_ k 

1  +  (k/x o  -  l)e~er 


□ 


(1.4) 


Since  e~01  0  as  t  -*•  oo,  it  follows  from  Equation  (1.4)  that  x{t)  ->  k  as  t  -»  oo. 

Also,  since  0  <  x0  <  k,  the  denominator  is  always  greater  than  1,  so  the  population 
size  is  less  than  k  for  all  t.  So  k  is  the  carrying  capacity  of  the  habitat. 

Usually  the  situation  being  modelled  is  such  that  x0  is  much  smaller  than  k ;  for 
example,  when  a  species  is  introduced  into  a  habitat,  the  initial  population  size  is 
small  compared  to  the  carrying  capacity  of  the  habitat.  For  this  situation, 

Figure  1.1  shows  the  graph  of  x(t),  the  population  size  at  time  t  for  the  logistic 
model;  the  graph  is  called  the  logistic  curve.  For  small  t ,  the  curve  is 
approximately  exponential;  for  larger  t,  the  rate  of  growth  levels  off  and  then 
decreases,  the  height  of  the  curve  approaching  the  value  k  as  t  -*  oo. 


Figure  1.1  The  logistic  curve 
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Question  1.2 

(i)  Show  that  when  x0>k>0,  the  solution  of  Equation  (1.3),  the  differential 
equation  for  the  logistic  model,  is  given  by  Expression  (1.4).  Sketch  the  graph 
of  x(t )  against  t  in  this  case. 

(ii)  Find  the  solution  of  Equation  (1.3)  when  x0  =  k.  □ 

The  solution  to  this  question  shows  that,  whatever  the  initial  size  of  the 
population,  after  a  sufficiently  long  time  the  population  size  will  be  approximately 
k.  If  x0  >  k,  the  population  size  will  decrease  monotonically  towards  the  value  k;  if 
0  <  x0  <  k,  the  population  size  will  increase  monotonically  towards  the  value  k; 
and  if  x0  =  k,  the  size  will  remain  at  that  value.  That  is,  x  =  k  is  a  stable 
equilibrium  point  of  the  model:  whatever  the  size  of  the  population  at  a  given  time, 
it  will  thereafter  approach  the  value  k.  This  means  that  if  the  population  size  were 
to  increase  above  the  value  k  (for  example,  due  to  introductions  from  outside), 
then  it  would  immediately  begin  to  decrease  towards  k.  Or  if  some  disaster  were  to 
cause  a  sudden  reduction  in  the  population  size,  then  the  population  would 
immediately  begin  to  recover,  its  size  once  again  approaching  the  value  k. 

For  the  models  discussed  in  this  section,  the  birth  and  death  rates  depend  only  on 
the  size  of  the  population;  it  is  assumed  that  the  birth  rate  is  independent  of  the 
ages  of  the  individuals  in  the  population,  and  that  each  individual  is  as  likely  to 
die  at  any  one  age  as  at  any  other.  However,  in  reality  birth  rates  are  different  for 
females  of  different  ages  and  death  rates  are  different  for  different  age  groups.  The 
models  discussed  in  the  remaining  sections  of  this  unit  take  this  into  account. 


2  Life  tables 


Bills  of  Mortality  for  London  were  first  published  early  in  the  sixteenth  century.  At 
first,  these  recorded  only  deaths  from  the  plague,  but  later  they  were  expanded  to 
include  christenings  as  well  as  plague  deaths  and,  towards  the  end  of  that  century, 
deaths  from  other  diseases  were  also  included.  But  it  was  only  in  the  seventeenth 
century  that  men  began  to  study  populations  from  a  scientific  point  of  view.  The 
Englishman  John  Graunt  (1620-74)  was  probably  the  first  to  do  so.  Despite  the 
unreliable  nature  of  the  data  in  the  Bills  of  Mortality,  Graunt  made  an  exhaustive 
study  of  the  information  contained  in  them.  In  1662  he  published  a  book  entitled 
Natural  and  Political  Observations  on  the  Bills  of  Mortality.  Among  the  topics  he 
discussed  was  the  sex-ratio:  he  noted  that  the  proportion  of  male  births  is  slightly 
greater  than  the  proportion  of  female  births,  an  observation  that  seems  to  have 
been  new  and  surprising  in  1662.  He  also  discussed  the  general  mortality  from 
different  kinds  of  diseases,  and  attempted  to  develop  a  mortality  table  of  the  type 
used  by  insurance  companies  today.  He  compiled  the  first  life  table:  in  his  book  he 
gave  the  proportions  of  people  dying  in  each  of  a  number  of  age  groups,  and  from 
this  deduced 

that  of ...  100  conceived  there  remains  alive  at  six  years  end  64 
At  Sixteen  years  end  40 


At  Twenty  six  25 

At  Tirty  six  16 

At  Fourty  six  10 

At  Fifty  six  6 

At  Sixty  six  3 

At  Seventy  six  1 

At  Eighty  0. 


This  simple  table  is  the  forerunner  of  the  more  detailed  life  tables  which  are  used 
today. 

In  Subsection  2.1,  the  idea  of  a  life  table  is  used  as  the  basis  for  the  development 
of  a  population  model  which  allows  for  different  death  rates  at  different  ages.  The 
nature  of  the  information  contained  in  a  life  table  is  discussed  briefly  in 
Subsection  2.2;  and  in  Subsection  2.3  an  example  is  given  to  show  how  a  life  table 
is  used  to  estimate  death  rates  for  individuals  of  various  ages. 
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2.1  The  life  table  function 


A  life  table  gives,  for  various  ages,  the  proportion  of  the  individuals  born  into  a 
population  who  are  still  alive  at  that  age.  For  example,  in  the  population  for 
which  John  Graunt  constructed  the  first  life  table,  only  40  out  of  every  100  born 
into  the  population  were  still  alive  at  age  16. 

Associated  with  any  life  table,  and  hence  with  the  population  it  describes,  is  a 
function  called  the  life  table  function;  in  a  stochastic  model  this  is  simply  a 
function  which  gives,  for  each  age  x,  the  probability  that  any  individual  will  be  still 
alive  at  age  x.  We  shall  suppose  that  the  lifetimes  of  the  individuals  in  some 
population  are  independent  and  identically  distributed,  and  that  the  length  of  life 
of  an  individual  born  into  the  population  is  a  positive  random  variable  X.  Then 
the  life  table  function,  <2(x),  is  defined  by 

Q(x)  =  P(X>x),  x  >  0. 

So  Q(x)  is  the  probability  that  the  length  of  life  of  an  individual  exceeds  x,  i.e.  that 
the  individual  survives  to  be  still  alive  at  age  x.  In  particular,  since  the  length  of 
life  of  an  individual  is  a  positive  random  variable,  that  is  P(X  >  0)  =  1,  we  have 

e(0)  =  1. 

A  quantity  of  importance  in  insurance  is  the  additional  length  of  time  that  an 
individual  aged  x  can  be  expected  to  live;  this  is  denoted  by  ex  and  for  x  >  0  it  is 
given  by  the  conditional  expectation 

E{X  -  x\X  >  x)  =  E(X\X  >  x)  -  x. 

Of  particular  interest  is  the  expectation  of  life  at  birth,  e0  =  E(X).  This  can  be 
calculated  directly  from  the  life  table  function,  as  you  are  now  asked  to  show. 

Question  2.1 

(i)  Express  the  life  table  function  Q(x)  in  terms  of  F(x),  the  cumulative 
distribution  function  of  the  random  variable  X. 

(ii)  Hence  show  that  the  expectation  of  life  at  birth  is  given  by  the  formula 

1*00 

e0  =  Q(x)dx.  □  (2.1) 

Jo 

Question  2.2  Use  Formula  (2.1)  to  find  the  expectation  of  life  at  birth  for  an 
individual  from  a  population  with  life  table  function  given  by 

Q(x)  =  (1  +  ocx)e~ax,  x  >  0, 
where  a  is  a  positive  constant.  □ 


(Note  that  for  the  life  table  function  of  Question  2.2,  Q(x )  is  approximately  zero  for 
large  x — a  property  that  would  be  required  when  modelling  a  human  population.) 

The  life  table  function  has  been  introduced  in  order  to  develop  a  model  which 
allows  for  different  death  rates  at  different  ages.  So  how  does  the  death  rate  vary 
with  age  in  a  population  with  life  table  function  Q(x)l  To  indicate  the  dependence 
of  the  death  rate  on  age  x,  a  function  h(x)  is  used  to  represent  the  death  rate  for 
individuals  aged  x.  This  means  that  the  probability  that  an  individual  aged  exactly 
x  will  die  in  the  next  short  time  interval  of  length  Sx  is  /t(x)<5x  +  o(<5x);  that  is, 

P(x  <Ar<x  +  ^x|AT>x)  =  h(x)  Sx  +  o(Sx).  (2.2) 

The  function  h(x)  is  called  the  age-specific  death  rate  of  the  population;  it  is  also 
sometimes  known  as  the  force  of  mortality. 


The  conditional  probability  on  the  left-hand  side  of  Equation  (2.2)  can  be  rewritten 
in  terms  of  the  life  table  function  Q(x): 


P(x  <2f<x  +  <5x|X>x)  = 


P(x  <  X  <  x  +  «5x) 
P(X  >  x) 


P(X  >x)-P{X>x  +  Sx) 
P(X  >  x) 


Q(x)  -  Q(x  +  Sx) 
Q(x) 


Q(x)  #  0. 


Note  that,  in  the  remaining  sections 
of  this  unit,  x  denotes  the  age  of  an 
individual,  whereas  in  Section  1  and 
in  previous  units,  x  denoted  the  size 
of  a  population. 

A  random  variable  is  positive  if  the 
only  values  it  can  take  are  positive. 
Thus,  we  are  considering  only  live 
births. 


When  you  study  Unit  13,  you  will 
see  that  h{x)  is  the  same  as  a  hazard 
function  in  the  theory  of  renewal 
models. 


Notice  that  h(x)  is  an  individual 
death  rate. 
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So 


Q(x)-Q(x  +  Sx ) 

- —q(x) - =  h(x)Sx  +  o(«5x), 

and  rearranging  this,  we  have 

ux\  = _ 1  Q(x  +  Sx)  -  Q(x)  o(<5x) 

Q(x)  Sx  +  Sx  ■ 

Letting  Sx  tend  to  zero,  we  obtain 

m*°-  (2.3) 

This  equation  expresses  the  age-specific  death  rate,  h(x),  in  terms  of  the  life  table 
function,  Q{x).  However,  sometimes  it  is  useful  to  express  the  life  table  function  in 
terms  of  the  age-specific  death  rate.  In  the  next  question  you  are  asked  to  derive 
an  expression  for  Q(x)  in  terms  of  h(x). 


Question  2.3  By  integrating  Equation  (2.3),  show  that 
Q(x)  =  exp(  —  f  h(u)du\  □ 


(2.4) 


Example  2.1 


In  the  more  elementary  models  for  population  growth  discussed  in  Units  7  and  8, 
such  as  the  pure  death  process,  the  immigration-death  process  and  the  simple 
birth-death  process,  each  individual  alive  at  time  t,  independently  of  all  other 
individuals,  has  a  probability  v  St  +  o(St)  of  dying  in  the  next  short  time  interval  of 
length  St.  No  account  is  taken  of  the  age  of  the  individual;  it  is  tacitly  assumed 
that  the  age-specific  death  rate  h(x)  has  the  constant  value  v.  Using  Equation  (2.4) 
to  find  the  life  table  function  Q(x)  corresponding  to  a  constant  age-specific  death 
rate  v,  we  obtain 


Q(x)  =  exp  I  - 


fx 


vdt  ]  =  e 


The  corresponding  c.d.f.  F{ x)  of  the  length  of  life  X  is 

m  =  i  -  Q(x) 


=  1  X  >  0, 

and  hence  the  corresponding  p.d.f.  f(x)  is 
f(x)  =  ve~vx,  x>  0. 

So  the  assumption  of  a  constant  death  rate  v  implies  that  the  length  of  life  of  an 
individual  member  of  the  population  has  an  exponential  distribution  with 
parameter  v,  thus  showing  that  the  only  lifetime  distribution  for  which  the  death 
rate  is  independent  of  age  is  the  exponential  distribution. 


The  mean  of  the  exponential  distribution  with  parameter  v  is  1/v,  so  when  the 
death  rate  is  v,  the  mean  length  of  life  of  an  individual  member  of  the  population 
is  1/v.  This  can  alternatively  be  obtained  directly  from  the  life  table  function  by 
using  Formula  (2.1)  for  the  expectation  of  life  at  birth.  This  gives 

Q(x)dx=  I  e~vxdx  =  -.  □ 

Jo  v 
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Question  2.4 

(i)  Find  the  life  table  function  corresponding  to  an  age-specific  death  rate  given 
by 

h(x)  =  fi2x,  x  >  0, 
where  fi  is  a  positive  constant. 

(ii)  Hence  find  the  p.d.f.  of  the  lifetime  distribution  of  an  individual.  Identify  this 
distribution,  and  hence  write  down  the  expectation  of  life  at  birth.  □ 


Question  2.5  Repeat  Question  2.4  for  a  population  in  which  the  age-specific  death 
rate  is  given  by 


h(x)  = 


1 


100 -x 
0 


0  <  x  <  100 
x  >  100.  □ 


Question  2.6  Find  the  life  table  function  and  the  age-specific  death  rate  of  a 
population  in  which  the  lifetime  X  of  an  individual  has  a  gamma  distribution  with 
parameters  (2,  a).  □ 


2.2  Life  tables 


Although  examples  such  as  the  ones  in  Subsection  2.1  are  useful  for  developing  an 
understanding  of  life  table  functions  and  age-specific  death  rates,  the  life  table 
functions  of  human  populations  do  not  have  forms  that  may  be  expressed  in  terms 
of  simple  mathematical  functions.  A  life  table  function  is  a  continuous  function 
but,  in  practice,  its  value  is  given  only  at  a  finite  sequence  of  ages;  for  example,  for 
human  populations  its  value  may  be  given  at  yearly  intervals  or  at  five-yearly 
intervals.  These  values  are  usually  presented  in  tabular  form;  that  is,  in  a  life  table. 

Clearly,  it  is  not  possible  to  tabulate  a  life  table  function  for  an  infinite  number  of 
ages  and,  in  any  case,  a  model  for  a  population  is  necessarily  based  on  a  finite 
amount  of  data.  So  values  are  estimated  and  tabulated  for  a  finite  set  of  ages,  and 
linear  interpolation  is  used  to  estimate  the  life  table  function  for  values  in  between. 
But  how  are  the  tabulated  values  of  a  life  table  function  obtained?  In  an  ideal 
world,  a  population  of  newborn  individuals  (all  born  at  approximately  the  same 
time)  would  be  observed  at  a  succession  of  times  and,  at  each  time  point  x,  the 
number  still  alive  would  be  recorded.  From  these  values  the  proportions  of  the 
population  who  have  survived  to  be  still  alive  aged  x  would  be  calculated  for  each 
x;  these  would  be  the  estimated  values  of  the  life  table  function.  A  life  table 
prepared  in  this  way  is  called  a  cohort  life  table. 


The  main  part  of  Graunt’s  life  table 
was  based  on  ten-yearly  intervals. 


In  Subsection  2.3  an  example  is 
given  to  show  how  a  life  table  can 
be  used  to  estimate  the  expectation 
of  life  at  birth  and  the  age-specific 
death  rate;  linear  interpolation  is 
used  implicitly  in  these  calculations. 


However,  in  the  real  world,  following  the  fate  of  a  cohort  is  difficult  for  human 
populations  (and,  indeed,  for  any  species  which  has  a  reasonably  long  lifespan  or 
which  cannot  be  reared  in  laboratory  conditions).  In  practice,  information  is 
collected  over  a  short  period  of  time  (as  is  the  case  for  census  data)  and 
observations  are  made  for  a  population  containing  individuals  of  all  ages.  A  life 
table  constructed  from  observing  an  all-age  population  is  called  a  current  life  table. 

For  the  present  purpose  of  studying  models  for  populations,  you  do  not  need  to 
know  the  details  of  how  life  tables  are  obtained,  or  how  the  life  table  function  can 
be  estimated  from  data  collected  by  observation  of  an  all-age  population.  Although 
the  life  tables  in  this  unit  have  not  been  obtained  directly  from  data,  they  are  all 
cohort  life  tables  and  they  will  be  referred  to  simply  as  life  tables.  You  do  not  need 
to  know  how  they  have  been  obtained,  but  you  should  be  aware  that  there  are 
difficulties  involved  in  estimating  the  life  table  function  for  a  human  population. 


2.3  Examples  of  life  tables 


At  any  given  time,  life  tables  vary  from  region  to  region;  and  within  a  region  they 
change  over  time  as  the  level  of  health  care  rises.  In  addition,  men  and  women 
have  substantially  different  life  tables.  Table  2.1  provides  an  example  of  a  life  table 
for  the  females  in  a  region;  it  is  a  model  table  for  the  females  in  a  human 
population.  It  will  be  used  to  illustrate  how  the  expectation  of  life  at  birth  and 
age-specific  death  rates  can  be  estimated  from  a  table.  The  first  two  columns 
constitute  the  life  table.  The  figures  in  the  third  column,  which  is  incomplete,  are 
called  central  mortality  rates;  they  are  calculated  from  the  information  in  the  first 
two  columns  and  are  frequently  included  in  a  life  table.  Later  in  this  subsection 
you  will  see  how  they  are  calculated  and  that  they  may  be  regarded  as 
approximations  to  values  of  the  age-specific  death  rate. 


Table  2.1  A  life  table  for  the  females  of  a  human  population 


Age  in  years, 

X 

/*  =  105  x  Q(x) 

Central  mortality  rate, 
mx 

0 

100000 

0.0415 

1 

95931 

0.0032 

5 

94693 

0.0010 

10 

94  212 

0.0008 

15 

93  838 

20 

93  244 

25 

92429 

30 

91468 

35 

90349 

40 

88  990 

45 

87  267 

50 

84931 

55 

81653 

60 

77061 

65 

70377 

70 

60737 

75 

47434 

80 

31368 

85 

16017 

90 

5  348 

95 

893 

100 

48 

105 

0 

This  table  was  estimated  from  data 
and  used  to  model  the  females  in  a 
population,  hence  it  is  a  model  table. 


This  table  is  taken  from  a  book  of 
demographic  tables:  A.  J.  Coale  and 
P.  Demeny,  Regional  Model  Life 
Tables  and  Stable  Populations 
(Academic  Press,  2nd  edn,  1983). 

The  calculations  carried  out  here  are 
of  a  simpler  kind  than  the  ones 
whose  results  are  listed  by  Coale 
and  Demeny. 


When  drawing  up  a  life  table,  in  order  to  avoid  the  occurrence  of  fractions  and 
decimals,  instead  of  tabulating  values  of  the  life  table  function  Q(x),  it  is  usual  to 
tabulate  values  of  100  x  Q(x),  for  example,  or  1000  x  Q(x)  or  100000  x  Q(x),  as 
appropriate.  In  Table  2.1,  values  of  100000  x  Q{x)  are  listed;  these  are  denoted  by 
lx,  so  lx  =  105  x  Q(x).  So  out  of  every  100000  females  born,  lx  is  the  expected 
number  who  are  still  alive  at  age  x. 

When  dealing  with  large  numbers  of  individuals,  random  fluctuations  may  be 
ignored,  and  a  deterministic  interpretation  of  the  life  table  may  be  adopted.  In  this 
case,  Q{x)  is  the  proportion  of  females  still  alive  at  age  x;  and  lx  is  the  number,  out 
of  every  100000  females  born,  who  are  still  alive  at  age  x.  From  now  on,  unless 
otherwise  stated,  a  deterministic  interpretation  will  be  adopted  in  this  unit. 

In  Table  2.1  the  values  of  lx  are  listed  at  five-yearly  intervals  with,  in  addition,  the 
value  for  x  =  1.  The  age  x  =  1  is  included  because  of  the  special  importance  of 
infant  mortality:  death  rates  in  the  early  months  of  life  are  much  higher  than  in 
the  following  years  of  childhood.  The  last  listed  age  is  105:  no  individual  survives 
to  age  105  years. 


The  expectation  of  life  at  birth 

The  life  table  may  be  used  to  calculate  an  approximate  value  of  the  expectation  of 
life  at  birth  which,  from  Formula  (2.1),  is  given  by  the  integral  of  the  life  table 
function  Q{x).  In  Figure  2.1  the  values  of  the  life  table  function,  Q(x)  =  lx  x  10“ 5, 
are  plotted  for  the  values  of  x  listed  in  Table  2.1.  The  plotted  points  are  linked  by 
straight  line  segments,  since  the  life  table  function  is  approximated  by  linear 
interpolation  for  intermediate  values  of  x. 
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Figure  2.1  The  life  table  function 


Figure  2.2  Approximating  the  area 


Now  the  integral  of  the  life  table  function  Q(x)  is  equal  to  the  area  under  the 
graph  of  Q(x),  so  the  expectation  of  life  at  birth,  e0,  may  be  approximated  by  the 
area  under  the  graph  in  Figure  2.1.  This  area  may  be  calculated  by  first  drawing 
vertical  lines  through  the  plotted  points,  thus  dividing  the  area  into  trapeziums,  as 
shown  in  Figure  2.2.  The  area  of  the  trapezium  between  x  and  x  +  k,  where  x  +  k 
is  the  next  age  listed  after  x,  is  equal  to  $h(lx  +  lx+k)  x  10" 5.  For  example,  the 
area  of  the  trapezium  from  x  =  85  to  x  +  k  =  90  is  \  x  5  x  (/85  +  /90)  x  10" 5;  this 
trapezium  is  illustrated  in  Figure  2.3. 

So  e0,  the  expectation  of  life  at  birth,  which  is  approximately  equal  to  the  total 
area  under  the  graph  in  Figure  2.1,  is  given  by 

eo  =  YJ^{lx  +  lx+k)  x  10" 5, 
x  L 

the  summation  being  over  all  listed  values  of  x.  For  x  =  0  and  x  =  1,  the  values  of 
k  are  1  and  4  respectively;  k  is  5  for  all  the  other  values  of  x.  So 

105  x  e0  =  -(/0  +  /j)  +  -(/!  +  /5)  +  -(/5  +  /10)  +  ^(fio  +  /15)  +  ...  +  |(/100  +  f i os) 

15  9  5 

=  2l°  +  2li  +  2^5  +  5  x  (^io  +  *15  +  •••  +  Jioo)  +  2*105 

=  6754266. 


Figure  2.3 


So  the  expectation  of  life  at  birth  for  females  in  this  population  is  approximately 
67.5  years. 


Question  2.7 

(i)  Plot  the  life  table  function  for  Graunt’s  life  table,  which  is  reproduced  in 
Table  2.2. 

(ii)  Use  the  life  table  to  calculate  an  approximate  value  of  the  expectation  of  life 
at  birth.  □ 


The  central  mortality  rate 

The  age-specific  death  rate  h(x)  is  given  by  Formula  (2.3): 


h{x)  = 


1  dQ 
Q(x)  dx ' 


To  use  this  formula,  the  life  table  function  is  required  in  the  form  of  a 
mathematical  function  that  can  be  differentiated.  The  life  table  function  as  drawn 
in  Figure  2. 1  is  not  differentiable  at  the  plotted  points,  where  two  straight-line 
segments  meet,  and  the  values  of  Q(x)  for  values  of  x  between  plotted  points  are 
only  approximations  obtained  by  linear  interpolation.  So  how  can  values  of  the 
age-specific  death  rate  be  obtained  when  the  life  table  function  is  given  only  for  the 
ages  listed  in  a  life  table? 


Table  2.2 


Age  in  years, 

X 

/,  =  100  x  Q(x) 

0 

100 

6 

64 

16 

40 

26 

25 

36 

16 

46 

10 

56 

6 

66 

3 

76 

1 

80 

0 
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One  approach  is,  for  each  listed  age  x,  to  calculate  a  value  mx  which  approximates 
h(x),  in  the  following  way.  First  the  life  table  function,  Q,  is  approximated  over  the 
range  x  to  x  +  k  by  the  average  of  its  values  at  x  and  x  +  k.  For  the  population  of 
Table  2.1,  this  is  equal  to 

(3(*)  +  Q(x  +  k)  _  lx  +  L+k 

2  “  2  x  105 ' 

Since,  in  Figure  2.1,  the  plotted  points  are  joined  by  straight-line  segments,  this 
value  is  equal  to  the  value  of  the  plotted  life  table  function  midway  between  x  and 

x  +  k;  that  is,  at  the  centre  of  the  age  range  x  to  x  +  k.  The  derivative  —  is 

dx 

approximated  by  the  gradient  of  the  line  segment  between  the  ages  x  and  x  +  k\ 
that  is  (for  Table  2.1),  by 

Q(x  +  k)-  Q(x)  ^lx+k-  lx 
k  k  x  105 ' 

So  mx,  the  approximate  value  for  the  age-specific  death  rate  h(x),  is  given  by 

2  x  105  lx+k  -  lx 

m,  =  —  - - - —  x 

lx  +  lx+k  fcxlO5’ 

that  is, 


m 


X 


2(1  x  ~  lx+k) 

k(lx  +  lx+k) 


(2.5) 


For  each  listed  age  x,  Formula  (2.5)  may  be  used  to  calculate  approximate  values 
mx  of  the  age-specific  death  rate;  mx  is  called  the  central  mortality  rate.  For 
example,  for  the  population  whose  life  table  is  given  in  Table  2.1, 


m0  = 


2(/0  ~  M 

Wo  +  W 


0.0415, 


2(1 1  ~  l5) 
4(/i  + 


=  0.0032, 


2(1 5  ~ 

Ws  +  ^lo) 


=  0.0010, 


m  io 

and  so  on. 


W10-/15) 
5(^o  +  hs) 


0.0008, 


The  word  central  reflects  the  fact 
that  the  approximation  is  based  on 
the  estimated  value  of  Q(x)  at  the 
centre  of  the  age  range  x  to  x  +  k. 


Question  2.8  Complete  the  column  of  values  of  the  central  mortality  rate,  mx,  for 
Table  2.1.  □ 


More  sophisticated  calculations  of  mortality  can  be  carried  out,  but  the  relatively 
crude  measure  defined  in  Formula  (2.5)  is  sufficient  to  show  how  death  rates  vary 
with  age.  The  values  of  the  central  mortality  rate,  mx,  are  plotted  against  x  (up  to 
70)  in  Figure  2.4. 

As  you  can  see,  after  the  relatively  high  mortality  of  infancy,  death  rates  are  at 
their  lowest  in  the  age  range  of  5  to  15  years.  Mortality  increases  slowly  during 
early  adulthood,  but  in  later  years  it  increases  ever  more  rapidly,  and  after  the  age 
of  70  it  is  large  enough  to  disappear  off  the  top  of  Figure  2.4. 

Question  2.9  Calculate  values  of  the  central  mortality  rate  for  Graunt’s  life  table 
(Table  2.2).  □ 


Figure  2.4  The  central  mortality 
rates  for  Table  2.1 
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Question  2.10 

(i)  Find  an  approximate  value  of  the  expectation  of  life  at  birth  for  a  population 
with  the  life  table  given  in  Table  2.3. 


Table  2.3  A  life  table  for  the 
females  of  a  human  population 


X 

1 X 

X 

K 

0 

100000 

45 

22  696 

1 

57180 

50 

20180 

5 

43  167 

55 

17318 

10 

40031 

60 

13  824 

15 

38403 

65 

9  756 

20 

36317 

70 

5  662 

25 

33  754 

75 

2464 

30 

30944 

80 

670 

35 

28  102 

85 

123 

40 

25  307 

90 

10 

95 

0 

(ii)  Calculate  values  of  the  central  mortality  rate  for  this  life  table. 

(iii)  Compare  your  results  for  this  life  table  with  those  for  the  life  table  in 

Table  2.1.  □ 


3  Stationary  populations 


In  this  section  and  the  next,  we  investigate  how  the  life  table  function  and  the 
growth  and  age-structure  of  a  population  are  related.  As  already  mentioned,  males 
and  females  have  different  life  tables;  in  fact,  women  tend  to  live  longer  than 
men — they  have  greater  expectations  of  life.  As  a  consequence,  the  distributions  of 
the  ages  of  the  males  and  the  females  in  a  population  are  different.  We  could  study 
the  age-distribution  of  either  males  or  females,  or  of  both  together  but,  since  the 
age-distribution  of  females  is  more  critical  for  the  birth  rate  and  hence  for  the  rate 
of  growth  of  a  population,  attention  is  restricted  to  the  size  and  age-distribution  of 
the  female  population.  Hence,  reference  to  a  life  table  will  imply  a  life  table  for 
females,  and  reference  to  birth  rates  will  imply  the  rates  at  which  daughters  are 
born  to  women. 

Incidentally,  it  is  generally  observed  that  the  sex-ratio  at  birth — that  is,  the  ratio 
of  male  births  to  female  births — is  greater  than  one,  typically  for  many 
populations  about  1.04  or  1.05.  Hence,  the  birth  rate  for  females  is  somewhat  less 
than  half  of  the  total  birth  rate. 

In  this  section  and  the  next,  we  consider  an  idealized  population  for  which  the  life 
table  function  does  not  change  over  time.  This  means  that  the  population  is 
assumed  to  be  living  over  a  long  period  of  time  in  conditions  which  undergo  little 
change  and  it  is  unaffected  by  any  major  catastrophe  such  as  a  famine,  a  plague  or 
a  war.  Furthermore,  it  is  assumed  that  the  population  is  closed;  that  is,  there  is  no 
immigration  or  emigration.  Again  this  is  an  idealization,  approximating  to  a 
situation  in  which  migration  plays  a  negligible  part  in  the  growth  of  the 
population. 


It  is  assumed  that  there  are  always 
enough  males  in  the  population  to 
maintain  the  birth  rate. 


Graunt  observed  that  this  was  so  in 
his  book  published  in  1662. 
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In  this  section  we  study  a  deterministic  model  for  a  population  which  is  constant 
in  size  and  whose  age-distribution  remains  constant  over  time;  that  is,  a 
population  that  does  not  change  in  size — births  are  balanced  by  deaths— and  in 
which  the  proportion  of  the  population  in  any  given  age  range  remains  fixed.  Such 
a  population  is  called  a  stationary  population.  In  Subsection  3.1,  we  investigate  the 
relationship  between  the  p.d.f.  of  the  age-distribution  of  a  stationary  population 
and  the  life  table  function  for  the  population.  And  in  Subsection  3.2,  various 
calculations  are  carried  out  for  stationary  populations  with  specific  life  tables. 


3.1  The  stationary  age-distribution 


In  a  stationary  population,  the  age-distribution  does  not  change  with  time  and  the 
size  of  the  population  remains  fixed,  so  births  must  occur  at  a  constant  rate.  We 
consider  a  stationary  population  of  constant  size  n  in  which  births  occur  at  a 
constant  rate  of  b  births  per  year.  The  per  capita  birth  rate  of  the  population  is 
b/n;  this  is  called  the  crude  birth  rate  of  the  population.  Of  course,  since  births  are 
balanced  by  deaths  in  a  stationary  population,  the  crude  death  rate  (that  is,  the 
per  capita  death  rate)  is  equal  to  the  crude  birth  rate. 


By  definition,  the  proportion  of  a  stationary  population  aged  between  xx  and  x2 
remains  constant  over  time  for  any  x,  and  x2  such  that  0  <  xx  <  x2,  so  the  age- 
distribution  of  the  population  may  be  represented  by  a  probability  density  function 
which  is  independent  of  time.  The  age-distribution  of  a  stationary  population  is 
called  the  stationary  age-distribution,  and  its  p.d.f.  is  denoted  by  g(x).  So  the 


proportion  of  the  population  aged  between  xt  and  x2  years  is  given  by  g(x)dx, 

Jxi 

and  the  total  number  of  individuals  in  the  population  aged  between  X!  and  x2  is 

rx  2 


equal  to  n 


g{x)dx. 


*i 


Clearly,  the  age-distribution  of  a  population  depends  on  the  life  table  function  Q{x) 
of  the  population;  for  example,  if  0(100)  =  0,  then  all  individuals  die  by  age  100, 
so  g(x)  =  0  for  x  >  100.  But  how  are  the  p.d.f.  g(x)  and  the  life  table  function  Q(x) 
related?  And  does  the  p.d.f.  also  depend  on  b,  the  rate  at  which  births  occur  in  the 
population,  or  on  the  crude  birth  rate  b/n. ?  To  answer  these  questions,  we  begin  by 
finding  two  expressions  for  the  number  of  individuals  in  the  population  at  time  t, 
who  are  aged  between  x  and  x  +  Sx;  for  convenience  we  call  this  number  N.  One 
of  these  expressions  involves  the  p.d.f.  g(x)  and  the  other  involves  the  life  table 
function  0(x). 


Firstly,  by  the  definition  of  the  p.d.f.,  N  is  given  by 

fx  +  dx 

N  =  kJ  g(y)dy 

so,  for  small  Sx,  we  have  N  ~  ng(x)Sx,  and  we  may  write 

N  =  n  g(x)  Sx  +  o(Sx).  (3.1) 

To  obtain  an  expression  for  N  involving  the  life  table  function  Q(x),  it  is  helpful  to 
use  a  type  of  diagram  known  as  a  Lexis  diagram;  an  example  is  shown  in 
Figure  3.1.  An  individual  aged  x  at  time  t  is  represented  on  the  diagram  by  a  point 
(x,  t).  The  collection  of  all  such  points  for  an  individual  from  her  birth  to  her  death 
is  her  life  line,  a  diagonal  line  starting  at  (0,  t0),  where  t0  is  the  time  of  her  birth, 
and  ending  at  (x^O,  where  tx  is  the  time  of  her  death  and  xx  =  tx  -  t0  is  her  age 
at  death.  (If  the  same  scale  is  used  on  each  axis,  the  life  line  is  at  45°  to  the  age 
axis.) 
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Time 


The  Lexis  diagram  in  Figure  3.2  will  be  used  to  obtain  an  expression  for  N 
involving  the  life  table  function  Q{x).  The  life  line  of  an  individual  who  is  aged 
between  x  and  x  +  dx  at  time  t  intersects  the  horizontal  line  at  time  t  between 
ages  x  and  x  +  dx;  Figure  3.2  shows  one  such  life  line.  An  individual  aged  x  at 
time  t  was  born  at  time  t  —  x,  and  an  individual  aged  x  +  dx  at  time  t  was  born  at 
time  t  —  x  —  dx;  the  life  lines  of  such  individuals  are  shown  by  the  dashed  lines  on 
the  diagram.  You  can  see  that  the  individuals  who  are  aged  between  x  and  x  +  dx 
at  time  t  are  those  who  were  born  in  the  time  interval  [£  —  x  —  dx,  t  —  x]  and  who 
have  survived  to  be  still  alive  at  time  t.  The  number  born  in  any  time  interval  of 
duration  <5x  is  b  dx,  and  the  proportion  of  those  born  in  the  interval 
[£  —  x  —  <5x,  t  —  x]  who  are  still  alive  at  time  t  is  approximately  Q{x)  so,  for  small 
dx,  we  have 

N  =  b  Q(x )  dx  +  o(<5x).  (3.2) 

From  Expressions  (3.1)  and  (3.2),  we  deduce  that 

ng(x)dx  =  bQ(x)dx  +  o(<5x). 

Rearranging  this  and  letting  dx  -+  0,  we  obtain 

g(x)  =  ^Q(x).  (3.3) 


So  the  p.d.f.  of  the  stationary  age-distribution  is  proportional  to  the  life  table 
function,  the  constant  of  proportionality  being  equal  to  the  crude  birth  rate  b/n. 
However,  the  crude  birth  rate  itself  depends  on  the  life  table  function — for 
example,  it  will  be  larger  for  a  population  in  which  most  individuals  die  young 
than  for  a  population  of  the  same  size  in  which  most  individuals  are  long-lived — 
so  how  is  it  related  to  Q(x)7  Since  g(x)  is  a  p.d.f.,  we  have 

•cO 

g(x)dx  =  1. 

Jo 


So  integrating  Equation  (3.3),  we  obtain 


Q(x)dx, 


and  using  Result  (2.1)  gives 


so 


1 


(3.4) 


That  is,  the  crude  birth  rate  is  equal  to  the  reciprocal  of  the  expectation  of  life  at 
birth.  Substituting  Result  (3.4)  in  Equation  (3.3),  we  find  that  the  p.d.f.  of  the 
stationary  age-distribution  is  equal  to  the  life  table  function  divided  by  the 
expectation  of  life  at  birth: 


g(x)  =  Q{x)/e0,  x  >  0. 


(3.5) 


An  immediate  consequence  of  Result  (3.5)  is  that,  for  a  stationary  population  with 
known  e0,  the  life  table  function  can  be  obtained  directly  from  the  p.d.f.  of  the  age- 
distribution  of  the  population.  So  if  it  is  reasonable  to  assume  that  a  given 
population  is  stationary  (or  at  least  approximately  stationary),  then  the  life  table 
function  can  be  estimated  very  simply  from  the  data  collected  on  a  population 
containing  individuals  of  all  ages;  it  is  not  necessary  to  observe  a  cohort  of 
individuals  over  a  long  period  of  time. 


Example  3.1 

Suppose  that  a  stationary  population  has  life  table  function  given  by 

<2(x)  =  (1  +  ax)e~ax,  x  >  0,  (3.6) 

where  a  =  0.05.  What  proportion  of  such  a  population  is  under  15  years  old,  or 
aged  65  and  over? 


t-x 

t  -  x  -  dx 

Figure  3.2 


Solution 


To  answer  this  we  must  first  find  the  age-distribution  of  the  population  using 
Result  (3.5).  The  expectation  of  life  at  birth  may  be  obtained  using  Formula  (2.1): 

r<x> 


e o  = 


Q(x)dx.  However,  from  Solution  2.6,  we  know  that  Q(x)  is  the  life  table 


Jo 

function  of  X  ~  T(2,a),  which  has  mean  2/a,  so  e0  =  2/a.  Therefore,  using 
Result  (3.5),  the  p.d.f.  of  the  stationary  age-distribution  is  given  by 


g{x)  =  |(1  +  xx)e  ax,  x  >  0, 


that  is,  since  a  =  0.05  =  1/20, 

?w=^(1+l))e""20’  xi0- 


The  proportion  of  the  population  under  15  years  old  is  given  by 


15 


^Wiix  =  40 


■15 


1  +  ^  )e  x,2°  dx. 


Using  integration  by  parts  to  evaluate  this  integral  gives  0.3505,  so  approximately 
35%  of  the  population  are  under  15  years  old.  Similarly,  the  proportion  of  the 
population  aged  65  and  over  is  given  by 


•oo 

•>65 


g(x)dx  =  - 


•oo 

65 


e~xl20dx  =  0.1018, 


that  is,  approximately  10%. 


A  statistic  of  social  and  economic  significance  for  any  population  is  the  ratio  of  the 
number  of  individuals  under  15  or  aged  65  and  over  to  the  number  aged  15  and 
over  but  under  65;  this  is  called  the  dependency  ratio.  For  the  population  of  this 
example,  the  dependency  ratio  is  equal  to 

(0.3505  +  0. 101 8)/(l  -  0.3505  -  0.1018)  =  0.83. 

So,  roughly  speaking,  in  this  case,  for  every  100  members  of  the  population  of 
working  age,  there  are  83  ‘dependent’  individuals.  □ 


Question  3.1  Suppose  that  a  stationary  population  has  life  table  function  given  by 
Formula  (3.6)  with  a  =  0.04.  Find  the  values  of 

(i)  the  expectation  of  life  at  birth; 

(ii)  the  proportion  of  the  population  under  15  years  old; 

(iii)  the  proportion  of  the  population  aged  65  and  over; 

(iv)  the  dependency  ratio.  □ 


Question  3.2  A  stationary  population  has  life  table  function  given  by 
Q(x)  =  1  -  x/100,  0  <  x  <  100. 

(i)  Find  the  crude  birth  rate  of  the  population. 

(ii)  Find  the  p.d.f.  of  the  stationary  age-distribution. 

(iii)  What  proportion  of  the  population  is  under  15  years  old? 

(iv)  What  proportion  of  the  population  is  aged  65  and  over? 

(v)  Find  the  dependency  ratio.  □ 


3.2  Stationary  populations  and  life  tables 

In  this  subsection  we  look  again  at  the  examples  of  life  tables  given  in  Subsection 
2.3,  and  in  addition  we  assume  that  each  is  the  life  table  of  a  stationary 
population.  In  Subsection  3.1  you  saw  that,  for  a  stationary  population  with  a 
given  life  table  function  Q(x),  the  crude  birth  rate  (and  hence  also  the  crude  death 
rate)  and  the  p.d.f.  of  the  age-distribution  can  be  calculated  directly  from  the  life 
table  function  using  Results  (3.4)  and  (3.5).  In  the  next  example  these  results  are 
used  to  calculate  approximate  values  for  the  crude  birth  rate  and  the  dependency 
ratio  for  a  stationary  population  with  the  life  table  in  Table  2.1. 


Page  10 
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Example  3.2 


By  Result  (3.4),  for  a  stationary  population  the  crude  birth  rate  b/n  is  equal  to 
l/<?0,  where  e0  is  the  expectation  of  life  at  birth.  In  Subsection  2.3,  for  Table  2.1, 
we  found  that  e0  ~  67.5  years,  so  the  crude  birth  rate  is  1/67.5  =  0.0148,  or 
approximately  15  per  thousand  population. 


The  expectation  of  life  at  birth  was  estimated  by  dividing  the  area  under  the  graph 
of  the  life  table  function  into  trapeziums  and  adding  their  areas.  A  similar 
approach  can  be  used  to  estimate  the  dependency  ratio. 


The  proportion  of  the  population  under  the 


•15 

Jo 


g(x)dx  =  — 
e0J 


•15 

0 


Q(x)  dx. 


age  of  15  is  given  by 


An  approximate  value  for 


•15 

Q(x)  dx  is  given  by  the  shaded  area  in  Figure  3.3. 
Jo 


So  we  have 


Figure  3.3 


f15 

105 !  ewrfx  a i(/„  + ),)  +  !(/,  + 15)  +  |(z5  +  ;10)  +  |(Z10  +  z15) 

1  5  %  5 

=  2l°  +  2*1  +  5*10  +  2? 15 
=  lJ&WfZ.  1 

So  the  proportion  of  the  population  under  the  age  of  15  is  14.689475/f?0  =  0.2175; 
or  approximately  3  j  <5J0 

Similarly,  the  proportion  of  the  population  aged  65  and  over  is 

/*i05  pl05 

g(x)dx  =  —  Q(x)dx, 

J 65  e0  J65 

pi  05 

and  Q(x )  dx  is  approximately  equal  to  the  shaded  area  in  Figure  3.4. 

J  65 


The  unrounded  value  of  e0  is  used 
here. 


Figure  3.4 
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So 

f105  55  5 

10s  Q(x)dx  ~  -(/65  +  /70)  +  x(^o  +  ^75)  +  ...  +  x(^ioo  +  ^105) 

65  Z  Z  Z 

=  2^65  +  5(/70  +  /75  +  ...  +  ll00)  +  2^105 

=  985  167.5. 

Therefore  the  proportion  of  the  population  aged  65  and  over  is 
9.851  675/e0  =  0.1459,  or  approximately  15%. 

The  proportion  of  the  population  aged  15  and  over  but  under  65  is 
1 — 0nH75~—  0rf459~=  0.636~6,~and  therefore  the  dependency  ratio  for  the 
population  is  (0^-7^H=^T459)/U.6366~^''0^7087-t)r  approximately  0.5$.  Roughly 
speaking,  this  may  be  interpreted  as  saying  that  for  every  100  members  of  the 
population  of  working  age,  there  are  .^‘dependent’  individuals — children  under 
15  and  the  elderly.  □ 

Question  3.3  For  a  stationary  population  with  the  life  table  given  in  Table  2.3,  Page  13 
obtain  approximate  values  for 

(i)  the  crude  birth  rate  and  the  crude  death  rate; 

(ii)  the  proportion  of  the  population  under  15  years  old; 

(iii)  the  proportion  of  the  population  aged  65  and  over; 

(iv)  the  dependency  ratio.  □ 

In  real  life,  the  crude  birth  rate  and  the  crude  death  rate,  commonly  expressed  as 
rates  per  thousand  population,  are  calculated  for  populations  including  males  as 
well  as  females.  Very  roughly  speaking,  including  the  males  in  the  population 
doubles  the  population  size  and  doubles  the  numbers  of  births  and  deaths  as 
compared  with  the  female  population  alone.  So,  for  a  stationary  population  whose 
females  have  the  life  table  of  Table  2.1  (for  example),  both  the  crude  birth  rate  and 
the  crude  death  rate  are  approximately  0.015,  or  15  per  thousand  population,  as 
calculated  in  Example  3.2  for  the  female  population. 


4  Stable  populations 

In  this  section  and  the  previous  one,  the  models  discussed  are  for  an  idealized 
population  for  which  the  life  table  function  (and  hence  the  age-specific  death  rate) 
does  not  change  over  time.  In  Section  3  we  studied  a  deterministic  model  for  a 
population  whose  age-distribution  remains  constant  over  time  and  which  is 
constant  in  size — a  stationary  population.  In  this  section  the  second  of  these 
additional  restrictions  is  relaxed,  so  the  population  considered  may  grow  in  size  or 
decline,  but  again  we  study  a  population  whose  age-distribution  remains  constant 
over  time.  How  the  size  of  the  population  changes  will  depend  on  the  age-specific 
death  rate  h(x)  or,  equivalently,  on  the  life  table  function  Q(x). 

The  deterministic  model  considered  in  this  section  is  for  a  population  in  which 
both  the  age-distribution  and  the  age-specific  death  rate  h(x)  (and  hence  <2(x)) 
remain  constant  over  time.  Such  a  population  is  called  a  stable  population — it  is 
important  to  remember  that  the  adjective  ‘stable’  does  not  necessarily  apply  to  the 
population’s  size.  Since  the  size  of  the  population  may  change  over  time,  it  is  a 
function  of  time,  t.  The  number  of  females  in  the  population  at  time  t  will  be 
denoted  by  n(t). 

As  in  earlier  sections  of  the  unit,  the  population  is  assumed  closed  to  immigration, 
so  new  members  can  join  only  by  births.  In  order  for  a  population  to  remain 
stable,  births  must  occur,  or  else  the  population  will  age  and  die  out.  Furthermore, 
these  births  must  occur  in  a  certain  pattern  over  time,  in  order  that  the  age- 
distribution  of  the  population  should  remain  fixed.  In  Section  3,  the  overall  birth 
rate  in  a  stationary  population  was  denoted  by  b.  In  the  case  of  a  stable 


This  comment  is  only  approximately 
valid  because,  as  has  already  been 
stated,  the  sex-ratio  at  birth  is  not 
quite  equal  to  one,  and  the  life 
tables  for  men  and  women  are 
somewhat  different. 
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population,  this  overall  birth  rate  will  in  general  be  a  function  b(t)  of  time  t;  that 
is,  the  total  number  of  births  which  take  place  in  the  small  time  interval  [£,  t  +  St] 
is  b{t)5t  +  o(St).  In  Subsection  4.1  we  investigate  the  way  in  which  the  birth  rate 
b{t)  and  the  population  size  n(t)  change  with  time  for  a  stable  population  with  a 
given  life  table  function. 

In  Subsection  4.2,  various  problems  are  discussed  for  a  stable  population.  For 
example,  what  form  is  taken  by  the  p.d.f.  of  the  age-distribution  of  a  stable 
population?  The  subsection  ends  with  a  discussion  of  the  relationship  between 
stable  populations,  the  exponential  model  and  stationary  populations. 

In  Subsection  4.3,  the  model  is  refined  and  further  investigated  by  taking  into 
account  more  explicitly  the  source  of  the  births  in  the  population.  For  most 
human  populations,  in  any  one  year  a  woman  aged  20  is  more  likely  to  give  birth 
than  is  a  woman  of  40.  Thus,  other  things  being  equal,  a  population  with  a 
predominance  of  20-year-olds  is  likely  to  have  more  births  than  is  a  population 
with  a  predominance  of  40-year-olds.  In  Subsection  4.3  a  deterministic  stable 
population  model  is  developed  which  deals  with  this  feature  of  human  populations. 


4.1  The  growth  of  a  stable  population 


In  this  subsection  we  begin  our  investigation  of  a  stable  population  by  finding  an 
equation  satisfied  by  the  overall  birth  rate  b(t).  This  is  done  in  a  similar  manner  to 
what  was  done  in  Subsection  3.1  for  stationary  populations;  we  find  two  different 
expressions  for  the  number  in  the  population  at  time  t  who  are  aged  between  x 
and  x  +  <5x.  To  do  this,  some  further  notation  is  introduced. 


Since  the  age-distribution  of  the  population  remains  constant  over  time,  its  p.d.f.  is 
denoted  by  g(x).  So  the  proportion  of  the  population  aged  between  x  and  x  +  <5x 
at  any  time  t  is  g(x)  <5x  +  o(<5x),  or  approximately  g(x)  <5x,  and  the  number  of 
females  in  the  population  at  time  t  who  are  aged  between  x  and  x  +  <5x  is 
n(t)g(x)Sx  +  o(<5x)  ~  n(t)g(x)Sx. 

A  second  expression  for  this  number  can  be  found  by  an  argument  similar  to  the 
one  used  in  Subsection  3.1.  From  the  Lexis  diagram  in  Figure  4.1  (which  is  the 
same  as  Figure  3.2),  you  can  see  that  the  individuals  who  are  aged  between  x  and 
x  +  <5x  at  time  t  are  those  who  were  born  in  the  time  interval  [£  -  x  -  <5x,  t  -  x] 
and  who  have  survived  to  be  still  alive  at  time  t.  The  number  of  births  in  the  short 
time  interval  [t  -  x  -  dx,t  -  x]  is  approximately  b{t  -  x)Sx,  and  of  these  a 
proportion  Q{x),  approximately,  survives  to  be  still  alive  at  time  t.  So  our  second 
expression  for  the  number  aged  between  x  and  x  +  dx  at  time  t  is  b(t  -  x)Q(x)Sx. 
Thus,  since  the  two  expressions  must  be  equal, 

n(£)g(x)<5x  =  b(t  -  x)Q(x)Sx, 
and  hence  we  have 

n(t)g(x)  =  b(t-x)Q(x).  (4-1) 

Since  Q(0)  =  1,  setting  x  =  0  in  Equation  (4.1)  gives 
n(t)g(0)  =  b(t), 

so  that 

n(t)  =  b(t)/g(  0). 

Then  substituting  for  n(t)  from  Equation  (4.2)  into  Equation  (4.1) 
throughout  by  g(x)  gives 

b(t)  _b(t-x)Q(x)_ 
m  g(x) 
so  for  x  >  0, 

b(t)  =  b(t  —  x)  c(x), 

where 


(4-2) 

and  dividing 


(4.3) 


From  now  on,  for  simplicity,  we 
shall  often  omit  the  terms  o(<5r)  and 
o(«5x). 


Figure  4.1 


c(x)  =  Q(x)g(0)/g(x),  (4.4) 

a  function  of  x  only  (i.e.  not  of  t).  Equation  (4.3)  relates  the  overall  birth  rate  b(t) 
at  time  t  to  the  corresponding  rate  b(t  -  x)  at  the  earlier  time  t  -  x.  Therefore,  if 
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for  instance  b{ 0)  is  known,  and  if  c(x)  is  known  for  all  non-negative  x,  then  putting 
t  =  x  in  Equation  (4.3)  gives 

b(x)  =  b( 0)  c(x),  x  >  0, 

and  so  b(t)  can  be  found  for  any  t  >  0.  In  fact,  it  is  possible  to  say  more. 

It  is  a  fairly  straightforward  exercise  to  show  that  any  function  b(t)  which  satisfies 
an  equation  such  as  Equation  (4.3)  must  be  an  exponential  function.  To  see  that 
this  is  so,  we  first  set  t  =  y  +  dy  and  x  =  dy  in  Equation  (4.3),  and  then  set  t  =  y 
and  x  =  0  in  Equation  (4.3);  these  steps  give 

b{y  +  8y)  =  b{y)c(8y) 

and 

b(y)  =  b(y)c(  0); 

so  subtracting  the  second  of  these  equations  from  the  first  and  dividing  by  5y,  we 
have 


b(y  +  dy)  -  b{y)  =  c(fty)  -  c(0) 
dy  dy 

Then,  letting  dy  tend  to  zero,  we  obtain 
b'(y)  =  6  b(y), 


where  6  =  c'( 0),  which  is  a  real  constant.  This  differential  equation  has  the  solution 
b(y)  =  b( 0)  e6y  for  y  >  0,  so  the  birth  rate  b(t)  is  given  by 


b(t)  =  b(0)eet,  t  >  0; 
and  hence,  from  Equation  (4.2),  we  have 


n(t)  = 


m  =  m 

g(  0)  g(  0)' 


(4.5) 


or 

n(t)  =  n{ 0)  ee‘,  t  >  0.  (4.6) 

Hence,  in  any  stable  population,  both  the  overall  birth  rate,  b(t),  and  the  size,  n(t), 
of  the  population  grow  (or  shrink,  depending  on  the  sign  of  6)  exponentially  with 
time.  The  parameter  0  is  called  the  intrinsic  growth  rate  of  the  population. 

If  the  life  table  function  Q(x)  and  the  p.d.f.  g(x)  of  the  age-distribution  of  a  stable 
population  are  known,  then  they  can  be  substituted  in  Equation  (4.4)  to  give  c(x), 
and  the  intrinsic  growth  rate  0  can  then  be  found  using  0  =  c'(0).  But  in  practice  it 
is  not  usual  to  have  available  values  of  g(x)  that  can  be  used  in  this  way.  The 
problem  is  normally  the  other  way— given  Q(x)  and  0,  to  find  g(x).  (In  practice, 
£>(x)  and  9  are  estimated  from  the  available  data,  and  g(x)  is  found  using  them.  In 
this  course  we  do  not  consider  such  estimation.) 

The  value  of  0  clearly  depends  in  some  way  on  the  balance  between  births  and 
deaths.  Later  in  this  section  we  shall  investigate  how  the  value  of  0  is  related  to 
other  features  of  the  population.  But  for  the  present,  this  problem  will  be  shelved. 
Instead,  the  age-distribution  of  a  stable  population  will  be  found. 


4.2  Stable  populations,  stationary  populations  and  the  exponential  model 

The  age-distribution  of  a  stable  population  is  called  the  stable  age-distribution.  In 
the  following  question  you  are  asked  to  find  an  expression  for  the  p.d.f.  g(x)  of  this 
distribution. 

Question  4.1 

(i)  Starting  from  Equation  (4.3),  show  that  c(x)  =  e6x. 

(ii)  Hence  show  that  the  p.d.f.  g(x)  of  the  stable  age-distribution  satisfies 

g(x)  =  g(0)e-exQ(x),  x  >  0.  □ 


(4.7) 


The  value  of  g(0)  may  be  obtained  by  using  the  fact  that,  since  g(x)  is  a  p.d.f., 

/'cO 

g(x)dx  =  1.  Integrating  both  sides  of  Equation  (4.7)  between  the  limits  0  and  oo 
Jo 

leads  to 


g(  0)  = 


e~exQ{x)dx 


so  the  stable  age-distribution  is  given  by 
e~9xQ(x) 


g(x)  = 


e~0xQ{x)dx 


(4-8) 


Stable  populations  and  the  exponential  model 

Since  the  size  of  a  stable  population  increases  (or  decreases)  exponentially  with  the 
intrinsic  growth  rate  9,  it  appears  that  this  population  model  may  be  a  special  case 
of  the  exponential  model  which  was  discussed  at  the  beginning  of  Section  1.  In 
fact,  this  is  so,  as  will  now  be  demonstrated. 

The  exponential  model  is  a  deterministic  model  corresponding  to  the  simple  birth- 
death  process.  In  the  model,  the  overall  birth  rate  and  death  rate  are  proportional 
to  the  size  of  the  population,  so  that  the  per  capita  birth  rate  /?  and  the  per  capita 
death  rate  v  do  not  change  over  time.  In  this  model  the  intrinsic  growth  rate  9  is 
related  to  the  per  capita  birth  and  death  rates  /?  and  v  by  the  equation  9  =  /?  -  v. 
Are  the  overall  birth  and  death  rates  in  a  stable  population  proportional  to  the 
size  of  the  population?  As  you  will  see,  this  is  indeed  the  case,  so  that  the 
exponential  growth  of  a  stable  population  occurs  essentially  because  such  a 
population  fits  the  exponential  model  of  Section  1. 

Question  4.2  Show  that,  in  a  stable  population,  the  overall  birth  rate  b(t)  satisfies 
the  equation  b(t)  =  p  n(t)  for  all  t  >  0,  where  p  is  a  constant.  Show  that 
P  =  g(  0).  □ 

Solution  4.2  shows  that  the  overall  birth  rate  in  a  stable  population  is  indeed 
proportional  to  the  size  of  the  population,  so  that  the  per  capita  (or  crude)  birth 
rate  ft  does  not  vary  over  time.  It  is  slightly  more  complicated  to  show  that  the 
same  is  true  for  death  rates.  Suppose  that  the  overall  rate  at  which  deaths  occur  at 
time  t  is  d(t) — that  is,  the  number  of  deaths  that  occur  in  the  short  time  interval 
[t,t  +  <5f]  is  d(t)5t  +  o(5t).  We  require  an  expression  for  the  overall  death  rate  d(t ); 
this  is  the  subject  of  the  next  question. 

Question  4.3 

(i)  By  obtaining  an  expression  which  does  not  include  d(t)  for  the  number  of 
deaths  in  the  short  time  interval  [£,  t  +  5t],  show  that  the  overall  death  rate 
d(t)  of  a  population  with  age-specific  death  rate  h(x)  is  given  by 

J%00 

n(t)g{x )  h(x)dx. 
o 

(ii)  Hence  show  that  the  overall  death  rate  d{t)  satisfies  the  equation  d(t)  =  v  n(t), 
where  v  is  a  constant.  □ 

This  shows  that  a  stable  population  has  a  per  capita  (or  crude)  death  rate  v  that 
does  not  vary  with  time.  Since  the  per  capita  rates  ft  and  v  do  not  vary  with  time, 
the  discussion  of  the  exponential  model  in  Section  1  shows  that  the  population  size 
must  increase  (or  decrease)  exponentially  at  intrinsic  rate  9  =  ft  —  v. 


Stable  populations  and  stationary  populations 


In  the  special  case  when  the  intrinsic  growth  rate  of  a  stable  population  is  zero, 
the  population  size  remains  constant  over  time,  and  the  population  has  all  the 
properties  which  define  a  stationary  population.  So  the  stationary  population 
discussed  in  Section  3  is  a  special  case  of  a  stable  population  in  which  the  intrinsic 
growth  rate  9  is  zero.  In  this  case,  the  integral  in  the  denominator  of 

J*  CO 

Q(x)dx,  which  is  equal  to  e0,  the  expectation  of  life 
o 

at  birth.  So,  in  this  special  case,  from  Equation  (4.8),  g(x)  =  Q(x)/e0;  this  is 
Equation  (3.5),  which  relates  the  p.d.f.  of  the  age-distribution  of  a  stationary 
population  to  the  life  table  function.  So  we  may  regard  the  p.d.f.  of  the  stable 
age-distribution,  given  by  Equation  (4.8)  as 


£(*)  = 


e~exQ(x) 


e~exQ(x)dx 


as  a  generalization  of  the  p.d.f.  of  the  age-distribution  of  a  stationary  population 
with  the  same  life  table  function.  In  Equation  (4.8),  the  life  table  function  is 
weighted  by  the  exponential  function  e~6x.  This  implies  that,  if  9  >  0,  so  that  the 
population  is  increasing,  then  the  proportion  of  the  population  in  the  youngest  age 
groups  is  greater  than  that  in  the  corresponding  stationary  population;  but  if 
0  <  0,  so  that  the  population  is  decreasing,  then  the  proportion  in  the  youngest 
age  groups  is  smaller  than  that  in  the  corresponding  stationary  population. 
Furthermore,  the  greater  the  value  of  6,  the  larger  is  the  proportion  of  the  stable 
population  in  the  youngest  age  groups.  For  example,  for  a  given  life  table  function 
Q(x),  the  shape  of  the  p.d.f.  might  be  as  shown  in  Figure  4.2  for  the  cases  9  >  0  (a 
growing  population),  9  =  0  (a  stationary  population)  and  9  <  0  (a  shrinking 
population). 


Figure  4.2  The  p.d.f.  of  the  stable  age-distribution 


In  the  discussion  in  Section  5,  the  shape  of  the  age-distribution  of  each  of  two  real 
human  populations  is  used  in  conjunction  with  the  theory  of  stable  populations  to 
gain  insight  into  the  way  the  populations  are  developing.  But  before  that,  we  shall 
investigate  the  theory  of  stable  populations  a  little  more  deeply. 

It  is  worth  pausing  at  this  point  to  take  stock  of  what  has  been  shown  in  this 
section  so  far.  A  stable  population  is  one  in  which  the  life  table  function  and  the 
age-distribution  do  not  vary  with  time.  It  has  been  shown  that  in  such  a 
population,  the  overall  size  of  the  population  n(t)  and  the  overall  birth  rate  b{t) 
both  change  exponentially  over  time  at  the  same  rate  9 ,  called  the  intrinsic  growth 
rate.  It  has  been  shown  that  the  p.d.f.  of  the  age-distribution  of  the  population, 
g(x),  must  satisfy  Equation  (4.8).  Finally,  it  has  been  shown  that  the  population  has 
per  capita  birth  rate  fi  and  death  rate  v  that  do  not  vary  with  time,  and  that  they 
satisfy  9  =  0  -  v.  But  methods  have  not  yet  been  introduced  to  enable  the  intrinsic 
growth  rate  6  to  be  found  for  a  stable  population  whose  life  table  function  Q(x)  is 
known.  In  the  next  subsection  you  will  see  how  to  get  round  this  problem. 


4.3  The  intrinsic  growth  rate  of  a  stable  population 

The  size  of  a  stable  population  may  increase,  decrease  or  remain  constant 
depending  on  whether  the  intrinsic  growth  rate  9  is  greater  than,  less  than  or  equal 
to  zero.  Whether  the  population  grows  or  declines  depends  on  the  birth  and  death 
rates  in  the  population.  Specifically  in  terms  of  mortality,  it  depends  on  the  age- 
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specific  death  rate,  or  equivalently  on  the  life  table  function,  which  we  are 
assuming  to  be  fixed  over  time.  In  terms  of  births,  we  shall  need  an  analogue  of 
the  age-specific  death  rate  which  takes  account  of  the  different  propensities  of 
women  of  different  ages  to  give  birth.  This  feature  of  populations  can  be  modelled 
using  the  age-specific  birth  rate,  A(x),  of  the  population,  which  is  defined  as  the  rate 
at  which  females  of  age  x  give  birth  to  females.  So,  for  instance,  the  number  of 
daughters  born  to  a  female  aged  x  in  the  short  time  interval  [t,  t  +  St ]  is 
/[\)St  +  o(St).  In  general,  the  age-specific  birth  rate  will  vary  with  time  as  well  as 
with  age  x.  However,  in  many  human  populations  it  turns  out  to  vary  rather 
slowly  with  time,  and  such  populations  can  sometimes  be  modelled  well  by  a 
model  in  which  2(x)  is  constant  over  time.  In  this  subsection,  then,  we  consider  a 
stable  population  model  in  which  the  age-distribution,  the  age-specific  death  rate 
and  the  age-specific  birth  rate  all  remain  fixed  over  time. 


So  the  value  of  9  depends  on  the  life  table  function  Q(x)  and  the  age-specific  birth 
rate  A(x).  As  a  first  step  towards  finding  a  relationship  between  9  and  these  two 
given  functions,  we  shall  obtain  an  equation  for  the  birth  rate  b(t)  involving  Q(x) 
and  A(x),  but  not  including  the  p.d.f.  g(x).  We  can  then  go  on  to  substitute  into  this 
equation  the  formula  for  b(t)  obtained  in  Equation  (4.5). 

You  have  already  seen  that  the  number  of  females  in  the  population  who  are  aged 
between  x  and  x  +  Sx  at  time  t  is  n(t)g(x)Sx.  In  the  short  time  interval  [t,t  +  <5f], 
these  females  give  birth  to  n(t)g(x)Sx  x  X(x)St  daughters,  so  the  total  number  of 
daughters  born  to  females  of  all  ages  in  the  interval  [t,  t  +  SQ  is 

^  n(t)g(x)^(x)dxjdt.  But,  by  definition,  the  number  of  female  births  in  the 

interval  [t,t  +  <5r]  is  b(t)St,  and  hence 

b(t)=  j  n(t)g(x)X(x)dx. 
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From  Equation  (4.1),  n(t)g(x)  =  b(t  -  x)Q(x),  and  so 

foo 

b(t  —  x)  Q(x)  2(x)  dx. 


(4.9) 


Given  the  life  table  function  Q(x)  and  the  age-specific  birth  rate  2(x),  Equation  (4.9) 
is  an  integral  equation  for  the  unknown  function  b(t).  It  has  already  been  shown 
(Equation  (4.5))  that  the  birth  rate  at  time  t  is  given  by  b(t)  =  b(O)e0t,  so  if  we 
substitute  for  b(t)  and  b(t  —  x)  in  this  integral  equation,  then  we  have 


b(O)e0t  = 


•oo 

b{0)  e0(t  ~  x)Q{x)  A(x)  dx ; 
Jo 


and  if  we  divide  each  side  by  b( 0)  e0t,  then  we  obtain  the  following  equation  for  the 
intrinsic  growth  rate  9: 

f*co 

e~0xQ(x)X(x)dx  =  1.  (4.10) 
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The  intrinsic  growth  rate  9  is  a  root  of  this  equation.  However,  this  equation  does 
not  provide  an  explicit  expression  for  9,  although  for  a  given  life  table  function 
and  age-specific  birth  rate,  iterative  methods  may  be  used  to  evaluate  solutions  of 
the  equation. 


Question  4.4  In  this  question  you  are  asked  to  investigate  the  general  kind  of 
process  involved  in  solving  Equation  (4.10)  for  9,  by  solving  it  in  a  simple  but  very 
unrealistic  case. 

(i)  Suppose  that  the  product  Q(x), l(x)  is  equal  to  1  for  x  between  0  and  2,  and  is 
equal  to  0  otherwise.  (This  is  not  meant  to  be  a  realistic  model  for  any  human 
population!)  Show  that  any  solution  of  Equation  (4.10)  must  satisfy  the 
equation 

1  -e~20  =  9.  (4.H) 

(ii)  Use  Newton-Raphson  iteration  with  a  starting  value  9  =  1  to  solve  Equation 
(4.11)  to  an  accuracy  of  three  places  of  decimals.  □ 


It  is  not  necessary  for  2(x)  to  be 
constant  over  time  for  the 
population  to  be  stable,  as  long  as 
the  overall  birth  rate  b(t)  varies  over 
time  in  the  appropriate  way.  Thus, 
the  model  considered  in  this 
subsection  is  a  special  case  of  the 
model  considered  earlier  in  the 
section. 


The  age-specific  birth  rate  2(x)  is 
zero  outside  the  fertile  age  range  for 
women,  so  the  range  of  integration 
could  in  practice  be  taken  as  being 
from,  say,  x  =  15  to  x  =  50. 
However,  we  shall  continue  to  write 
the  limits  of  integration  as  x  =  0  to 

X  =  OO. 
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If  realistic  functions  Q(x)  and  A(x)  are  used,  the  process  of  solving  Equation  (4.10) 
for  8  is  generally  more  involved  than  what  you  did  in  Question  4.4,  but  there  are 
iterative  and  approximate  methods  which  can  be  used.  (We  shall  not  go  into  the 
details  of  such  calculations  in  this  course.) 

Is  it  possible  for  Equation  (4.10)  to  have  more  than  one  solution?  To  answer  this 

J%00 

e~wxQ(x)X(x)dx.  From 
o 

Equation  (4.10),  the  intrinsic  growth  rate  8  is  a  value  of  w  that  is  a  root  of  the 
equation  /(w)  =  1.  By  considering  properties  of  the  functions  e~wx,  Q(x)  and  2(x) 
(including  the  property  that  A(x)  is  non-zero  only  on  a  finite  interval— the  fertile 
age  range),  it  can  be  shown  that  the  general  shape  of  the  graph  of  /(w)  against  w  is 
that  shown  in  Figure  4.3:  that  is,  /(w)  is  a  continuous  monotonic  decreasing 
function  of  w  such  that  /(w)  co  as  w  ->■  -  oo  and  /(w)  ->  0  as  w  -►  oo.  Thus  the 
graph  of  I(w)  must  cross  the  horizontal  line  corresponding  to  the  value  1  once  and 
only  once.  Hence,  the  equation  /(w)  =  1  has  a  unique  real  root,  and  this  root  must 
therefore  be  the  intrinsic  growth  rate  8  of  the  stable  population. 

For  a  stable  population,  the  birth  rate  b(t),  and  hence  also  the  population  size  n(t), 
either  increases  or  decreases  exponentially.  But  is  the  intrinsic  growth  rate  8 
positive,  negative  or  zero?  Does  the  population  size  increase  or  decrease  with  time, 
or  remain  constant,  as  is  the  case  for  a  stationary  population?  This  can  be 
answered  by  looking  at  the  graph  of  /(w)  against  w.  The  graph  crosses  the  vertical 
axis  at  R0,  where  R0  is  given  by 

R0  =  m=\  Q{x)X{x)dx.  (4.12) 
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The  value  of  R0  is  determined  by  the  life  table  function  <2(x)  and  the  age-specific 
birth  rate  A(x)  and,  depending  on  the  particular  Q(x)  and  A(x),  it  may  be  greater 
than,  equal  to  or  less  than  one.  In  Figure  4.4(a),  the  graph  crosses  the  vertical  axis 
at  R0>  1,  and  in  this  case  8  >  0.  In  Figure  4.4(b),  it  crosses  the  vertical  axis  at 
R0  <  1,  so  that  in  this  case  8  <  0. 


Figure  4.4 


In  the  special  case  R0  =  1,  we  have  0  =  0,  and  the  population  size  remains 
constant,  so  that  the  population  is  stationary.  So,  depending  on  the  given  life  table 
function  and  age-specific  birth  rate,  the  population  size  can  either  increase 
exponentially,  remain  constant  or  decrease  exponentially.  Which  of  these  occurs 
can  be  determined  by  using  Equation  (4.12)  to  calculate  the  value  of  the  constant 
R0,  which  is  called  the  net  reproduction  rate  of  the  population.  The  reason  for  this 
name  should  become  apparent  when  you  do  the  next  question. 


Question  4.5  In  a  stochastic  interpretation  of  the  age-specific  birth  rate,  the 
expected  number  of  female  offspring  born  to  a  female  of  age  x  in  the  next  short 
time  interval  of  length  <5x  is  A(x)  Sx  +  o(Sx). 

(i)  Using  this  interpretation,  show  that  the  expected  total  number  of  female 
offspring  born  to  a  female  during  her  lifetime  is  given  by 

co 

Ro=  Q(x)A{x)dx. 

Jo 

(ii)  Explain  on  intuitive  grounds  why  8  >  0  if  R0  >  1  and  8  <  0  if  R0  <  1.  □ 


If  it  were  possible,  we  would  have  to 
investigate  which  of  the  solutions 
gave  the  intrinsic  growth  rate  of  the 
population. 


/(w) 


w 


Figure  4.3  The  graph  of  the 
function  /(w) 
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In  this  subsection  we  have  shown  that,  in  a  stable  population  whose  age- 
distribution,  age-specific  birth  rate  and  age-specific  death  rate  are  all  fixed  with 
time,  the  intrinsic  growth  rate  must  be  the  unique  real  root  9  of  Equation  (4.10). 
Therefore,  given  the  life  table  function  Q(x )  and  the  age-specific  birth  rate  A(x)  for 
such  a  population,  one  can  find  the  intrinsic  growth  rate,  and  by  substituting  it  in 
Equation  (4.8)  the  p.d.f.  g{x)  of  the  stable  age-distribution  can  be  found.  The 
practical  importance  of  such  a  procedure  stems  mainly  from  the  following  result, 
whose  proof  is  beyond  the  scope  of  the  course. 

Suppose  a  population  develops  over  time  in  such  a  way  that  its  life  table  function 
Q(x)  and  its  age-specific  birth  rate  2(x)  remain  fixed  over  time.  Then,  whatever  the 
initial  age-distribution  of  the  population,  its  age-distribution  will  approach  the 
stable  age-distribution  given  by  Equation  (4.8)  with  the  value  of  9  given  as  the  root 
of  Equation  (4.10).  In  the  long  run,  then,  any  population  whose  age-specific  birth 
and  death  rates  do  not  change  will  become  approximately  stable,  and  the  theory  in 
this  section  can  be  used  to  investigate  its  age-structure  and  growth  rate.  In 
practice,  many  human  populations  do  have  age-specific  birth  and  death  rates  that 
change  fairly  slowly  with  time,  so  that  their  development  over  time  can  be 
modelled  reasonably  accurately  by  treating  them  as  stable.  However,  this  is  not  the 
case  for  all  populations,  as  you  will  see  in  the  next  section. 


5  Population  pyramids 

In  the  previous  three  sections  of  this  unit,  models  have  been  developed  for  the 
growth  of  a  human  population.  In  this  section  two  examples  of  age-distributions 
are  discussed  qualitatively.  These  are  taken  from  the  published  results  of  the  1970 
census  of  the  USSR. 

The  USSR  comprises  fifteen  union  republics,  and  Figures  5.1  and  5.2  illustrate  the 
age-distributions  of  two  of  them — the  Uzbek  republic  and  the  Latvian  republic — 
at  the  time  of  the  census.  The  age-distributions  are  illustrated  by  diagrams  called 
population  pyramids.  These  are  histograms  which  plot  population  size  by  age 
group  as  a  percentage  of  the  total  population  size,  with  males  on  the  left-hand  side 
and  females  on  the  right-hand  side.  In  the  published  census  data,  ages  are  grouped 
in  five-year  intervals  up  to  the  age  of  sixty  and  in  ten-year  intervals  thereafter;  the 
same  grouping  is  used  in  Figures  5.1  and  5.2.  The  percentage  scale  along  the 
horizontal  axis  is  for  the  five-year  intervals. 


Figure  5.1  Population  pyramid  for  the  Figure  5.2  Population  pyramid 

Uzbek  republic  in  1970  for  the  Latvian  republic  in  1970 


Since  a  population  pyramid  is  a 
histogram,  the  area  of  each  bar  is 
proportional  to  the  percentage  of 
the  population  in  that  age  group. 
Thus,  the  length  of  a  bar  for  a  ten- 
year  interval  (in  Figures  5.1  and  5.2) 
is  half  the  percentage  in  that  age 
group. 
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A  fundamental  difference  between  the  shapes  of  the  population  pyramids  for  the 
Central  Asian  Uzbek  republic  in  Figure  5.1  and  for  the  European  Latvian  republic 
in  Figure  5.2  is  immediately  apparent.  The  pyramid  of  Figure  5.1  has  a  broad 
base,  showing  that  the  younger  age  groups  are  very  strongly  represented,  and 
hence  indicating  that  the  population  has  been  growing  rapidly.  On  the  other  hand, 
the  pyramid  of  Figure  5.2  has  a  much  narrower  base,  showing  that  the  younger 
age  groups  are  somewhat  less  represented  than  they  would  be  in  a  stationary 
population,  and  hence  at  first  sight  suggesting  that  the  population  may  even  have 
been  contracting  slightly. 

These  comments  have  to  be  qualified  by  the  observation  that  the  recent  history  of 
both  republics  has  been  such  that  the  conditions  assumed  in  the  development  of 
the  stable  theory  have  not  been  approached  even  approximately.  Levels  of  fertility 
have  changed  over  time,  the  populations  have  been  affected  significantly  by 
migration,  and  the  effects  of  the  Second  World  War  and  its  aftermath  are  still 
clearly  visible  in  the  population  pyramids.  The  detailed  study  of  the  age-structure 
of  the  populations  would  require  considerable  knowledge  of  the  history  of  the 
USSR  and  lies  in  the  domain  of  the  historical  demographer  rather  than  the 
mathematical  statistician.  Nevertheless,  some  of  the  effects  of  the  war  are  easy  to 
see.  For  instance,  in  Figure  5.1  for  the  Uzbek  republic,  the  relatively  small 
numbers  in  the  age  groups  from  20  to  29  years  old  are  a  consequence  of  the 
decrease  in  birth  rates  during  the  war  and  the  years  immediately  following.  In 
Figure  5.2  for  the  Latvian  republic,  this  feature  is  masked  by  the  effects  of 
immigration  on  a  large  scale.  During  the  nineteen-forties  the  Latvian  republic 
suffered  enormous  population  losses  due  to  war  deaths,  deportations  and 
emigration,  balanced  to  a  large  extent  by  post-war  immigration;  this  immigration, 
mainly  of  young  adults,  has  strongly  influenced  population  growth.  War  losses 
reveal  themselves  in  the  population  pyramids  of  Figures  5.1  and  5.2  through  the 
under-representation  of  the  male  population  in  the  age  groups  over  45  and  the 
especially  small  numbers  of  both  sexes  in  the  50-54  age  group. 

These  examples  show  up  the  limitations  of  stable  population  theory  and  may  lead 
us  to  question  its  relevance.  However,  the  mathematical  theory  does  help  to  give 
us  a  general  understanding  of  the  age-structure  of  populations  and  provides  us 
with  a  yardstick  against  which  to  measure  observed  population  pyramids.  In  some 
cases,  populations  will  have  a  more  or  less  stable  age-structure;  and,  where  this  is 
not  the  case,  the  difference  between  the  observed  structure  and  what  is  predicted 
by  the  stable  population  theory  may  provide  the  demographer  with  a  focus  for 
further  investigation. 


Objectives 

After  studying  this  unit  you  should  be  able  to: 

describe  the  exponential  model  and  the  logistic  model  for  the  growth  of  a 
population; 

derive  and  use  the  equations  relating  the  life  table  function  of  a  population  and  the 
age-specific  death  rate;  and  find  the  expectation  of  life  at  birth  for  a  population 
with  a  given  life  table  function; 

derive  and  use  the  equation  relating  the  life  table  function,  the  p.d.f.  of  the  age- 
distribution  and  the  expectation  of  life  at  birth  of  a  stationary  population; 

use  the  information  in  a  life  table  to  obtain  approximate  values  for  the  expectation 
of  life  at  birth  and  the  central  mortality  rates;  and,  for  a  stationary  population, 
obtain  approximate  values  for  the  crude  birth  rate  and  the  dependency  ratio; 

show  that  the  population  size  and  the  overall  birth  rate  in  a  stable  population 
change  exponentially  over  time; 

find  expressions  for  the  p.d.f.  of  the  age-distribution  and  the  per  capita  birth  and 
death  rates  of  a  stable  population  with  intrinsic  growth  rate  0; 

discuss  the  relationships  between  stable  populations,  stationary  populations  and 
the  exponential  model  for  a  population; 

obtain  the  integral  equation  satisfied  by  the  intrinsic  growth  rate  of  a  stable 
population  which  has  age-specific  birth  rates  constant  over  time. 
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Appendix:  Solutions  to  questions 

Section  1 


l.l(i)  Separating  the  variables  in  Equation  (1.3),  we  obtain 
k 


But 


x(k  —  x) 
k 


dx  =  \9dt. 


1  1 


x(k  —  x)  k  —  x  x 
so  we  have 


+  -  I  dx  =  I  Odt, 


v/c  —  x  x, 
giving  (since  k  >  x  >  0) 

—  log  {k  —  x)  +  log  x  =  6t  +  constant 
or 


k  -  x 

log  I  — - —  )  =  —  Qt  +  constant. 


(ii)  Using  the  initial  condition  x(0)  =  x0  gives 
k-x0 


constant  =  log 


x0 


so 


log 


k  —  x 


=  log 


k  —  x, 


X  J  -\  X0 
Talcing  exponentials,  we  have 


-Qt. 


*-1  =  (A_ne-., 

x  \x0 


which  leads  to 

x(t)  = 


1  +  {k/x0  -l)e~ 


t>  0. 


1*2(0  The  solution  of  Equation  (1.3)  when  x0  >  k  >  0 
differs  from  that  in  Solution  1.1  only  at  one  point.  When 
x  >  k,  integrating  l/(k  -  x)  leads  to  -log(x  -  k)  since  the 
domain  of  the  logarithm  function  is  the  positive  real 
numbers;  that  is,  logy  is  defined  only  for  y  >  0.  This  leads  to 
t  ,  n 

—  0t  +  constant,  x  >  k. 


log 


The  initial  condition  x(0)  =  x0  leads  to  Expression  (1.4)  as 
before.  (From  this  expression  it  follows  that  x(t)  >  k  for  all  t, 
and  so  there  is  no  need  to  consider  the  case  x  <  k.)  The 
graph  of  x(t)  against  t  is  sketched  below. 


dx  . 


(ii)  If  x0  =  k,  then  —  is  zero  at  time  0,  so  the  population 


size  will  not  change:  x(f)  =  k  for  all  t  >  0. 


Section  2 


2.1(i)  By  definition,  Q(x)  =  P(X  >  x),  so 
Q(x)  =  1  -  P(X  <  x)  =  1  -  F(x). 

(ii)  The  expectation  of  life  at  birth  is  e0  =  E(X).  Using 
Formula  (i^l)  from  Unit  1  for  the  mean  of  a  positive  random 
variable,  we  have 


e0=  \  (1  -  F(x))dx 

) 

<2(x)  dx,  using  part  (i). 


2.2  Using  Formula  (2.1),  we  have 

Pco 

e0  =  (1  +  ax)e~axdx 


=  [(1  +  ax)e-«/(-a)]y  + 


cdx, 


using  integration  by  parts, 
=  1/a  +  1/a  =  2/a. 


2.3  Integrating  Equation  (2.3)  (using  a  dummy  variable  u), 
we  have 


h{u)du  =  )0-<k)%du 


So 


=  [-logQ(u)]g,  since  Q(u)  >  0  in  (0,x), 
=  -logQ(x)  +  log  0(0) 

=  —  log  Q(x),  since  0(0)  =  1. 


Iog0(x)=  - 


h{u)  du. 


and,  taking  exponentials,  we  have 
0(x)  =  exp^-  f  h(u)dii). 


2.4(i)  Using  Formula  (2.4),  we  have 
0(x)  =  exp(^-  f  p2udu 


=  e  P2*2!2,  x  >  0. 

(ii)  The  c.d.f.  of  the  lifetime  X  is 
F(x)  =  1  -  0(x)  =  1  -  e-P2*2'2, 
so  the  p.d.f.  of  X  is 
fix)  =  F(x) 

=  p2xe~fiY12,  x  >  0. 

This  is  the  p.d.f.  of  a  Rayleigh  distribution  with  parameter  (1, 

and  this  has  mean  -  yJn/2.  So  the  expectation  of  life  at  birth 

.  1  / —  p 

is-y/n/2. 
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2.5(i)  Using  Formula  (2.4),  we  have,  for  0  <  x  <  100, 

e(x)  =  exp(-|Tj^du) 

=  exp  ([log  (100  -  m)]5) 

—  exp  (log  (1  —  x/100)) 

=  1  -  x/ 1 00, 

and  Q(x)  =  0  for  x  >  100. 

(ii)  The  c.d.f.  of  X  is 

F(x)  =  x/100,  0  <  x  <  100, 
and  hence  the  p.d.f.  is 

f{x)  =  1/100,  0  <  x  <  100. 

So  X  is  uniformly  distributed,  X  ~  U{ 0, 100),  and  has  mean 
50.  The  expectation  of  life  at  birth  is  50. 


2.6  Since  X  ~  T(2,  a),  its  p.d.f.  is 
/(x)  =  a2xe~ax,  x  >  0. 

The  c.d.f.  of  X  is 


F(x)  =  a 2ue  du 
Jo 

=  [a2ue““7(— a)]o  ~ 


u2e~auf (-ot)du. 


Jo 


using  integration  by  parts, 


=  -axe  ax  +  ae  au  du 
Jo 

=  —axe~ax  +  [  —  e_““]o 


=  —axe  —  e 


+  1. 


So  the  life  table  function  is 
Q(x)  =  1  -  F{x) 

=  (1  +  ax)e~ax,  x  >  0. 

(This  is  the  life  table  function  of  Question  2.2.) 
Using  Formula  (2.3),  the  age-specific  death  rate  is 

ea~ 


h{x)  =  - 


1  +  ax 

2 


x  (ae  ax  -  a(  1  +  ax)e  1X) 


ax 


1  +  ax’ 


x  >  0. 


2.7(i)  The  life  table  function  is  shown  below. 


(ii)  An  approximate  value  of  the  expectation  of  life  at  birth 
is  given  by 
k 

e0  =  \ ]j(lx  +  lx+k)  x  10-2, 

X  *■ 

SO 

102  X  e0  =  -(/0  +  lb)  +  —  (/6  +  he)  +  +  ^26)  +  ••  • 

10  4 

+  y(^66  +  he)  +  2^76  +  ^8o) 

=  3/0  -I-  8/6  +  10(/16  +  ...  +  l66)  +  ll16  +  2/80 
=  1819. 

So  for  the  population  on  which  Graunt  based  his  life  table, 
the  expectation  of  life  at  birth  was  approximately  18.2  years. 

2.8  The  complete  column  of  values  of  the  central  mortality 
rate,  mx,  is  given  in  the  table  below. 


X 

mx 

X 

mx 

0 

0.0415 

50 

0.0079 

1 

0.0032 

55 

0.0116 

5 

0.0010 

60 

0.0181 

10 

0.0008 

65 

0.0294 

15 

0.0013 

70 

0.0492 

20 

0.0018 

75 

0.0816 

25 

0.0021 

80 

0.1296 

30 

0.0025 

85 

0.1997 

35 

0.0030 

90 

0.2855 

40 

0.0039 

95 

0.3592 

45 

0.0054 

100 

0.4000 

2.9  The  values  of  mx  are  given  below. 


X 

0 

6 

16 

26 

36 

46 

mx 

0.073 

0.046 

0.046 

0.044 

0.046 

0.050 

56 

66 

76 

0.067 

0.100 

0.500 

For  example, 

m0  =  2(100  -  64)/[6(  1 00  +  64)]  =  0.073, 
m6  =  2(64  -  40)/[  10(64  +  40)]  =  0.046 
and 

m16  =  2(1  -  0)/[4(  1  +  0)]  =  0.500. 

Notice  that  the  lack  of  detailed  information  on  deaths  during 
infancy  means  that  it  is  not  possible  to  detect  the  high  rate 
of  infant  mortality  which  also  occurred  in  Graunt’s  time. 

2.10(i)  It  is  not  necessary  to  draw  the  graph  of  the  life 
table  function  in  order  to  calculate  the  areas  of  the 
trapeziums  under  the  graph.  The  areas  can  be  calculated 
directly  from  the  values  of  lx  given  in  the  life  table.  The 
approximate  value  of  e0  is  given  by 
10s  x  c0 

=  2^0  +  h)  +  +  h)  +  |(/ 5  +  ho)  +  •••  +^90  +  hs) 

15  9  5 

=  2^°  +  2!1  +  2^5  +  ^,0  +  •••  +  ho)  +  2^95 

=  2015  006.5, 

so  the  expectation  of  life  at  birth  is  approximately  20  years. 
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(ii)  The  central  mortality  rates  are  given  in  the  table 
below. 


X 

mx 

X 

mx 

0 

0.5449 

45 

0.0235 

1 

0.0698 

50 

0.0305 

5 

0.0151 

55 

0.0449 

10 

0.0083 

60 

0.0690 

15 

0.0112 

65 

0.1062 

20 

0.0146 

70 

0.1574 

25 

0.0174 

75 

0.2290 

30 

0.0193 

80 

0.2759 

35 

0.0209 

85 

0.3398 

40 

0.0218 

90 

0.4000 

(iii)  The  life  table  values  for  Table  2.3  are  lower  and  the 
mortality  rates  are  higher  than  those  for  Table  2.1.  The 
individuals  in  this  population  tend  to  die  much  younger.  As 
a  consequence,  the  expectation  of  life  at  birth  is  only  about 
20  years  compared  with  68  years  for  the  population  of 
Table  2.1. 


Section  3 


(ii)  The  p.d.f.  of  the  stationary  age-distribution  is 
g(x)  =  Q(x)/e0,  so 


g(x)  = 


0  <  x  <  100. 


l~mdx 


(iii)  The  proportion  of  the  population  under  15  years  old 
is  given  by 

g(x)dx  =  ^ 

~  50 
=  0.2775, 
or  approximately  28%. 

(iv)  The  proportion  of  the  population  aged  65  and  over  is 
•100  ^ 


x2  T 

200 Jo 


g{x)dx  = 


50 

=  0.1225, 

or  approximately  12%. 


T 

200 ]6 


oo 

65 


(v)  The  dependency  ratio  is 

(0.2775  +  0. 1 225)/(  1  -  0.2775  -  0.1225)  =  0.6667, 
or  approximately  0.67. 


3.1(i)  As  in  Example  3.1,  the  expectation  of  life  at  birth  is 
e0  =  2/a,  so  when  a  =  0.04,  e0  =  50  years. 


(ii)  The  proportion  of  the  population  aged  under  15  years 
is  given  by 


g(x)dx  =  — 
eo  Jo 


Q(x)dx 


(1  +  <xx)e~ttXdx 


=  f  Tc(l  +ax)e‘“7(-<x)]o5  + 
using  integration  by  parts, 
=  i[-(2  +  «)e-“]J!; 


so,  when  a  =  0.04,  the  proportion  aged  under  15  years  is 
1  -  0.7135  =  0.2865, 
that  is,  approximately  29%. 


(iii)  Similarly,  the  proportion  of  the  population  aged  65 
and  over  is  given  by 

•co  ] 

g(x)dx  =  2  +  ctx)e  ax]£,; 
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so,  when  a  =  0.04,  this  proportion  is  0.1708,  that  is, 
approximately  17%. 


(iv)  The  dependency  ratio  is 

(0.2865  +  0.1708)/(1  -  0.2865  -  0.1708)  =  0.8426, 
that  is,  approximately  0.84. 


3.2(i)  The  crude  birth  rate  is  b/n  =  \/e0,  e0  being  the 
expectation  of  life  at  birth.  Now 


So  the  crude  birth  rate  is  1/50  =  0.02,  or  20  per  thousand  of 
the  population.  (Alternatively,  Solution  2.5  may  be  used  to 
determine  e0.) 


3.3(i)  From  Solution  2.10(i),  the  expectation  of  life  at  birth 

is  20.150065,  or  approximately  20  years,  so  the  crude  birth 

rate  and  the  crude  death  rate  are  both  equal  to 

1/20.15  =  0.04963,  or  approximately  50  per  thousand  of  the 

population. 

(ii)  The  proportion  of  the  population  aged  under  1 5  years 
is 

f  g{x)dx  =  —  f  Q(x)dx, 

Jo  eo  Jo 

where  e0  ~  20.15  years.  Dividing  the  area  under  the  graph  of 
Q(x)  into  trapeziums  between  x  =  0  and  x  =  15,  we  obtain 

il  5 

Q(x)dx  ~  -l0  +  -/j  +  -J5  +  5/10  +  -/15 
=  683  364. 

Hence  the  proportion  aged  under  15  years  is  (using 
unrounded  values)  6.833  64/20.150065  =  0.3391;  that  is, 
approximately  34%. 

(iii)  Similarly,  the  proportion  aged  65  and  over  is 
g{x)dx  =  —\  Q(x)dx, 

^65  e0  J 65 

and 

f°°  5 

105  Q(x)dx  ^  -l65  +  5 l-j0  +  5/75  +  5/80  +  5 /85  +  5/90 

^65  1 

=  69035. 

So  the  proportion  of  the  population  aged  65  and  over  is 
0.69035/20.150065  =  0.0343;  that  is,  approximately  3%. 

(iv)  The  dependency  ratio  is 

(0.3391  +  0.0343)/(l  -  0.3391  -  0.0343)  =  0.5959, 
or  approximately  0.60. 
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Section  4 


4.1(i)  From  Equation  (4.3),  c(x)  is  given  by 
c(x)  =  b(t)/b(t  -  x). 

Using  Equation  (4.5)  to  substitute  for  b(t)  and  b(t  -  x)  gives 
c(x)  =  b(0)e9t/(b(0)ee<,~x)) 

=  e6x. 

(ii)  From  Equation  (4.4),  the  p.d.f.  g(x)  is  given  by 
g(x)  =  Q(x)g(0)/c(x). 

So,  using  part  (i),  g(x)  satisfies 
g(x)  =  g(0)e~exQ(x),  x  >  0, 
as  required. 

4.2  There  are  (at  least)  two  ways  of  doing  this.  First,  from 
Equations  (4.5)  and  (4.6), 

so  that  b{t)  =  Pn(t)  with  =  b(0)/n(0).  Then,  from 
Equation  (4.2)  with  t  =  0,  h(0)/n(0)  =  g(0),  so  that  /?  =  g(0), 
as  required. 

Alternatively  (and  more  easily,  if  you  happen  to  spot  the 
right  equation  to  begin  with),  Equation  (4.2)  can  be  used 
directly.  It  states  explicitly  that  b{t)  =  g(0)n(t). 


4-3(1)  From  the  definition  of  the  overall  death  rate  d{t),  the 
number  of  deaths  that  occur  in  the  short  time  interval 
[£,  t  +  St ]  is  approximately  d(t)  St. 


An  expression  not  involving  d(t)  for  this  number  of  deaths  is 
rather  more  complicated  to  obtain.  At  time  t  there  are 
n(f)g(x)<5x  individuals  aged  between  x  and  x  +  Sx,  and  a 
proportion  h{x)St  of  these  will  die  in  [i,  t  +  <5f].  So  the 
number  aged  between  x  and  x  4-  Sx  at  time  t  who  die  in 
[£,  t  +  <5r]  is 

n(£)g(x)<5x  x  h(x)St  =  n(t)g(x)h(x)SxSt. 

Hence  the  total  number  of  individuals  of  all  ages  who  die  in 

/  r®  \ 


[t,  t  -f  <5t]  is 


n(t)g(x)h(x)dx 


St.  Therefore 


d{t)St  = 
and  so 


n(t)g(x)h(x)dx  )<5t 


d{t)=  n(t)g(x)h(x)dx. 

Jo 

(ii)  From  the  last  line,  since  n(f)  does  not  depend  on  x. 


d(t)  =  n(t ) 


g{x)h{x)dx. 


Jo 


Thus,  d(t)  =  vn(t),  where  v  =  j  g{x)h(x)dx  (which  does  not 

Jo 

depend  on  t )  is  a  constant. 


4.4(i)  In  this  case  the  integral  on  the  left  of  Equation  (4.10) 

becomes  e~9xdx,  which  is  equal  to  (1  -  e~20)/6  for  0^0. 
Jo 

Therefore  Equation  (4.10)  says  that  (1  -  e~2e)/6  =  1,  so  that 
\-e  29  =  9,  as  required.  (When  0  =  0,  the  integral  on  the 
left  of  Equation  (4.10)  is  equal  to  2.) 


(n)  To  solve  Equation  (4.1 1)  by  Newton-Raphson 
iteration,  it  must  be  written  in  the  form  /(0)  =  0.  This  can  be 
done  by  writing  it  as 
1  -  e~2e  -0  =  0, 

with  /(0)  =  1  —  e~29  —  0.  Therefore 
/'(0)  =  2e~2e  —  1, 

so  the  Newton-Raphson  iteration  formula  is 

0  =e  l~e~29'-dn 
n+1  *  2e~2e"  —  1  ' 

With  starting  value  0!  =  1,  this  formula  gives  successively 
02  =  O.81444, 

03  =  0.79702, 

04  =  0.79681, 

so  that  the  root  is  0  =  0.797  to  three  decimal  places. 


(You  may  have  noticed  that  0  =  0  is  also  a  solution  of 
Equation  (4.11).  However,  it  is  not  a  solution  of 
Equation  (4.10),  since  the  left-hand  side  of  Equation  (4.10)  is 
2  when  0  =  0.) 


4.5(1)  Consider  a  female  who  has  just  been  born. 
Conditional  upon  the  event  that  she  survives  to  age  x,  the 
expected  number  of  female  offspring  born  to  her  between 
ages  x  and  x  +  Sx  is  A(x)<5x  +  o(Sx).  The  probability  that  she 
survives  to  age  x  is  Q(x).  Hence,  unconditionally,  the 
expected  number  of  female  offspring  born  to  her  between 
ages  x  and  x  +  Sx  is  Q(x)  (A(x)  <5x  +  o(<5x)),  which  can  be 
written  as  Q{x)X{x)Sx  +  o(<5x).  Thus,  the  expected  total 
number  of  offspring  born  to  her  during  her  lifetime  is 


Q(x)X(x)dx  =  R0. 


That  is,  the  expected  number  of  female  offspring  born  to  a 
female  during  her  lifetime  is  equal  to  the  net  reproduction 
rate  R0. 


(ii)  If  R0>  1,  then,  on  average,  each  female  has  more  than 
one  daughter,  and  the  population  and  number  of  births  must 
increase  with  time  in  the  long  run.  Hence,  0  >  0. 

On  the  other  hand,  if  R0  <  1  then,  on  average,  each  female 
has  less  than  one  daughter,  and  the  population  and  number 
of  births  must  decrease  with  time  in  the  long  run.  Hence 
0<  0. 
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